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1.0  SCOPE 


SOFTWARE  DESIGN  DOCUMENT  FOR  THE 
NAVY  STANDARD  SURF  MODEL  VERSION  3.2 


1.1  Identification 

This  Software  Design  Document  (SDD)  is  written  for  the  Navy  Standard  Surf  Model,  version 
3.2,  or  SURF  3.2.  SURF  3.2  will  be  called  SURF  hereafter.  SURF  provides  users  with  an  automated 
method  for  determining  surf  conditions  and  related  environmental  parameters.  The  model  produces 
a  standard  surf  forecast  such  as  breaker  height,  longshore  current,  and  a  modified  surf  index  (MSI) 
number,  which  are  Navy  requirements  for  littoral  operations  and  amphibious  landings  (see  Joint  Surf 
Manual).  The  first  operational  Navy  surf  forecasting  computer  model  was  developed  for  the  Fleet  in 
1988  (Earle,  1988)  to  supplement  the  manual  and  visual  techniques  developed  in  the  1950’s.  The 
manual  procedures  are  subjective  and  do  not  adequately  consider  shallow  water  effects  such  as  wave 
shoaling  and  refraction.  SURF  is  a  modem  numerical  model  and  has  been  validated  by  field  and 
laboratory  data  (Hsu  et  al.  2000  and  2002). 

1.2  Document  Overview 

This  OAML  SDD  describes  the  design,  structure,  and  scientific  aspects  of  the  Computer 
Software  Configuration  Item  (CSCI)  titled  SURF.  This  document  provides  a  detailed  summary  of  all 
Computer  Software  Units  (CSU)  or  subroutines,  input  file  formats,  output  file  formats,  and  user- 
specified  options.  The  SDD  is  divided  into  three  sections:  the  Preliminary  Design,  the  Architectural 
Design,  and  the  Detailed  Design. 

The  Preliminary  Design  section  describes  the  scientific  aspects  of  SURF  including  a  brief 
description  of  the  mathematical  formulation  and  theory  behind  the  model.  The  Architectural  Design 
section  outlines  the  structural  design  of  SURF  with  a  graphical  representation  of  the  CSU  calling 
sequence.  The  Detailed  Design  section  identifies  and  summarizes  the  operation  of  each  CSU 
including  detailed  listings  of  input  variables,  output  variables,  local  variables,  calling  routines,  and 
called  routines  and/or  called  functions. 

2.0  REFERENCED  DOCUMENTS 
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Earle,  M.D.,  Surf  Forecasting  Software  Improvements,  MEC  Systems  Corp.  contract  report: 

N00014-91-C-601 1  to  NRL  (  Formerly  Naval  Oceanographic  and  Atmospheric  Research 
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Earle,  M.D.,  Surf  Forecasting  Software  User’s  Manual,  NRL  ( Formerly  Naval  Oceanographic 
and  Atmospheric  Research  Laboratory,  NORDA)  Tech.  Note  352,  ASW  Oceanography 
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3.0 


PRELIMINARY  DESIGN 


3.1  CSCI  Overview 

SURF  is  a  parametric  one-dimensional  model  based  largely  on  the  work  of  Thornton  and 
Guza  (1983,  1986).  Thornton  and  Guza  developed  several  models  for  random  wave  processes 
including  a  wave  height  transformation  model  and  a  longshore  current  model.  These  models  contain 
both  numerical  and  analytical  solutions,  which  provide  cross-shore  distributions  of  various 
parameters  such  as  wave  height,  longshore  current  velocity,  and  wavelength.  However,  because 
SURF  is  one-dimensional,  certain  approximations  are  made:  (1)  straight  and  parallel  bottom 
contours,  (2)  depth-uniform  currents,  (3)  wave  heights  are  Rayleigh  distributed,  (4)  linear  wave 
theory  is  applicable,  and  (5)  directional  wave  spectra  are  narrow-banded  in  frequency  and  direction. 

The  model  is  designed  to  operate  in  a  variety  of  modes  to  provide  both  military  and  civilian 
users  with  local  surf  and  current  forecasts.  SURF  requires  three  pieces  of  information  to  perform 
calculations:  (1)  a  depth  profile,  (2)  a  directional  wave  spectrum,  and  (3)  wave  refraction 
information.  Each  of  these  required  data  sources  can  be  accessed  externally  or  generated  internally. 
This  design  allows  for  maximum  flexibility  when  using  SURF  to  generate  local  forecasts  where 
input  data  may  or  may  not  be  available.  The  following  subsections  outline  the  scientific  principles 
behind  SURF  and  the  inherent  fundamental  hydrodynamic  calculations  contained  in  the  model. 

3.1.1  Wave  and  Roller  Energy  Models 

As  waves  approach  the  coast,  the  frictional  effect  of  the  sea  floor  on  the  organized  orbital 
motion  of  water  particles  within  a  wave  causes  the  wave  to  break  or  spill.  The  flows  of  spilling 
breakers  can  be  separated  into  two  layers,  an  upper  layer  of  turbulent  energy,  which  rides  over  a 
lower  layer  of  energy  that  maintains  an  organized  oscillatory  wave  motion.  The  region  of  turbulent 
water  above  the  wave  is  termed  a  surface  roller.  The  original  idea  of  such  a  two-layer  system  was 
introduced  by  Longuet-Higgins  and  Turner  (1974)  (see  also  Svendsen  (1984  a,b)).  SURF 
incorporates  the  model  of  Lippman  et  al  (1995,  1996),  which  produces  results  consistent  with 
measurements  from  both  a  planar  and  a  barred  beach.  The  energy  associated  with  each  region  of 
interest  is  utilized  to  shoal  the  incoming  waves  and  drive  the  longshore  current.  The  energy  per  unit 
surface  area  in  a  wave  is  calculated  as: 


1  2 
Ew  —  0  P8  H  mu 
o 

where  p  is  water  density  and  g  is  the  acceleration  due  to  gravity,  ff*  is  the  root-mean-square  wave 
height.  The  energy  per  unit  area  associated  with  a  roller  is  given  as: 

Er  8^  h tan a 

where  c  is  the  phase  speed  of  the  wave,/is  the  zero  crossing  frequency,  Hb  is  the  height  of  the  wave 
at  breaking,  h  is  water  depth,  and  a  is  the  angle  the  roller  makes  with  the  body  of  the  wave.  A 
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default  value  of  5  degrees  is  used  for  the  roller  angle. 
3.1.2  Energy  Dissipation  in  the  Surf  Zone 


As  a  wave  propagates  across  the  surf  zone,  its  energy  is  dissipated  due  to  bottom  friction 
wave  breaking,  turbulence,  and  wave-current  interaction.  A  formulation  of  this  energy  dissipation  is 
given  by  the  energy  flux  equation: 


d(  EwC,,  cos  6  ) 
dx 


=  <£b> 


where  Ew  is  the  wave  energy,  cg  is  the  wave  group  velocity  and  6  is  the  wave  direction  relative  to 
shore  normal  (x  positive  offshore).  The  Right  Hand  Side  (RHS)  of  the  above,  equation,  «*>,  is  the 

ensemble  averaged  dissipation  function.  Thornton  and  Guza  (1983)  modeled  this  dissipation 
function  as: 


<  e»  >  =  7  Pgf |-J  H3  Pb(H)dH 
4  n 

where  B  is  an  empirical  coefficient,  and  pb(H)  is  the  probability  distribution  for  breaking  waves 
described  by: 


pb(H)  =  W(H)p(H) 

where  p(H)  is  a  Rayleigh  Distribution  of  wave  heights  and  W(H)  is  a  weighting  function  resulting  in 
a  weighted  Rayleigh  distribution.  Several  weighting  functions  W(H)  have  been  constructed  by 

various  authors,  the  weighting  function  applied  in  SURF  developed  by  Thornton  and  Guza  (1986)  is 
given  as: 


where  y  is  an  empirical  factor  determined  from  field  data  to  be  0.42,  h  is  the  water  depth  and  H is  the 

wave  height.  If  wave  roller  energy  is  considered  in  the  model,  the  modified  energy  flux  equation  is 
given  as: 

d(  EwCgCOSQ  )  d(  ErCCOS0  ) 

dx  dx  -  <e> 

and  the  dissipation  becomes  a  function  of  the  roller  term. 
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<£r>  =  4 PSf^cosa  \ H3 pb(H)dH 
4  h 

The  above  equation  is  solved  using  a  numerical  forward  stepping  and  convergence  scheme  to 
determine  wave  and  roller  energy  along  with  Hrms  values  at  each  point. 

3.1.3  Longshore  Current  Calculations 

When  waves  enter  the  surf  zone  at  an  angle,  the  shore-parallel  component  of  momentum 
inherent  to  wave  motion  drives  a  current  along  the  shore.  This  longshore  current  can  be  a  significant 
force  inside  the  surf  zone.  Calculation  of  the  current  velocity  is  based  on  radiation  stress  theory  (see 
Longuet-Higgins,  1970a,  1970b).  A  general  form  of  the  longshore  momentum  equation  is: 

where  p  is  the  water  density,  h  is  the  water  depth,  and  V  is  the  longshore  current.  The  first  term  on 
the  left  hand  side  is  the  radiation  stress  in  the  along  shore  direction  exerted  by  waves  on  the  water 
given  by: 


Ty  =  <€b> 


sin# 

c 


where  £b  is  the  dissipation  function  defined  in  the  previous  section,  c  is  wave  phase  speed,  and  0  is 
the  angle  of  wave  approach  with  respect  to  x  The  second  term  is  the  horizontal  mixing.  The 
horizontal  eddy  viscosity  p  is  modeled  after  Battjes  (1975). 

in  which  M  is  an  empirical  constant  equal  to  2.  The  third  term  is  the  mean  stress  due  to  bottom 
friction  given  by: 


T  y=pcfuV 

where  c/  is  the  bottom  friction  coefficient,  u  is  the  magnitude  of  the  near-bottom  horizontal  wave 
oibital  velocity,  and  V  is  the  longshore  current.  Linear  wave  theory  defines  the  near-bottom  wave- 
induced  orbital  velocity  as: 

-  EH 

U  T  sinh  (kh) 

where  H  is  the  wave  height,  T  is  the  wave  period  and  k  is  the  wave  number,  which  can  be  calculated 
using  the  dispersion  relation: 
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<J2  =  gktanh(kh) 


where  cxis  the  radian  wave  frequency  and  g  is  gravity.  The  longshore  current  equation  is  solved 
using  a  finite  difference  approach  after  wave  heights,  water  depths,  and  wave  dissipation  values  are 
calculated  at  each  cross-shore  grid  point  in  the  surf  zone. 

A  major  improvement  to  the  longshore  current  calculation  is  included  in  SURF.  Hsu  et  al. 
(2000)  showed  that  using  a  variable  bottom  friction  coefficient  formulation  in  the  longshore  current 
model  provides  a  significant  improvement  in  longshore  current  velocities.  The  depth  dependent 
bottom  friction  coefficient  function  is  defined  as 


f/W= 


0.0035 


0.0033 


'x>' 

h— 

2 


h(x) 


> 


where  x  is  the  offshore  distance,  h  is  the  local  water  depth,  and  Xb  is  the  distance  from  the  shoreline 
to  the  location  where  ten  percent  of  the  waves  are  breaking. 

3.1.4  Directional  Energy  Spectra 


SURF  allows  users  to  generate  surf  forecasts  using  two  different  directional  wave  energy 
spectrum  types.  The  user  can  choose  from  an  internally  generated  wave  spectrum  or  an  external 
wave  spectrum.  If  the  internally  generated  spectrum  is  selected,  a  modified  Pierson-Moskowitz 
(1964)  spectrum  is  calculated  based  on  sea  and  swell  conditions  defined  in  the  surf  model  input  file. 

3.1.5  Differences  Between  SURF  3.1  and  SURF  3.2 


There  are  many  code  improvements  in  SURF  3.2.  In  general,  the  prolog  of  each  routine  in 
the  software  was  changed  to  reflect  individual  changes  and  in  many  routines,  comments  were  added, 
clarified,  deleted  or  corrected,  where  necessary.  The  major  changes  in  the  software  are  as  follows: 

1.  The  lookup  table  in  routine  MDSRF2  was  corrected  to  the  values  in  the  Joint  Surf  Manual. 
There  had  been  seven  typographical  errors,  which  produced  inconsequential  errors  in  the 
modified  surf  index.  In  the  process  of  correcting  these  errors,  the  routine  was  streamlined  and 
documented  better  than  it  had  been. 

2.  At  the  request  of  Systems  Integration  Division,  Naval  Oceanographic  Office,  the  input  variable 
dstart,  i.e.  the  starting  depth  of  the  surf  computation,  was  eliminated  from  the  SURF  input  file. 
This  variable  is  now  automatically  calculated  in  the  model  from  information  in  the  depth 
profile,  the  depth  of  the  input  waves  and  the  refraction  option.  The  order  of  input  parameters 
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was  re-arranged  for  clarity. 

3.  The  output  file  of  SURF  had  traditionally  produced  a  profile,  or  listing,  of  wavelength  through 
the  surf  zone  but  never  a  listing  of  wave  direction;  although,  directional  variation  throughout 
the  surf  zone  is  more  operationally  relevant  than  variation  in  wavelength.  Therefore,  in  the 
detailed  output,  wavelength  has  been  replaced  by  wave  direction.  This  change  required  adding 
another  routine  (lin4)  and  making  minor  changes  to  several  other  routines.  Extra  lines  were 
added  in  the  summary  part  to  list  the  average  wavelength  and  the  depth  of  wave  input. 

4.  The  routines  related  to  equilibrium  profile,  starting  point,  percent  of  breaking  and  wave 
refraction  option  were  improved. 

5.  Several  variables,  constants  and  arrays,  which  are  passed  through  several  unnecessary  levels  of 
computation,  were  either  eliminated  or  appropriately  moved. 

6.  Superfluous  routines  associated  with  experimental  versions  of  SURF  were  eliminated.  These 
routines  allowed  an  alternative,  but  unnecessarily  complex,  calculation  of  longshore  current. 
Routines  b_head  and  b_detail  were  merged  into  the  summary  routine. 

7.  The  angle  used  for  initializing  the  wave  transformation  through  the  surf  zone  was  replaced  from 
the  direction  of  the  vertically  averaged  momentum  flux  to  the  energy-weighted,  mean  wave 
direction.  The  change  was  made  because,  by  definition,  the  former  angle  limits  wave  directions 
to  within  45  degrees  of  shore  normal,  which  was  found  to  be  an  unrealistic  constraint.  In 
changing  the  definition  of  the  initializing  wave  angle,  it  was  necessary  to  increase  the  bottom 
friction  coefficient  from  0.0030  to  0.0035  after  re-calibrating  the  model  using  field  data. 

8.  The  cross-shore  distance  over  which  bottom  slope  is  calculated  for  subsequent  estimates  of 
breaker  type  was  increased  from  approximately  6  feet  to  a  one-quarter  the  wavelength.  The 
maximum  distance  is  limited  to  100  ft  and  the  minimum  distance  is  set  to  10  ft.  A  routine 
(get_slope)  was  added  for  this  calculation.  This  change  has  effectively  reduced  the  unusually 
high  incidence  of  estimates  of  surging  breakers. 

9.  The  interpolation  scheme  in  the  routine  grd_frc  was  completely  re-written  to  reflect  a  bilinear 
interpolation  scheme,  which  is  simpler  and  more  efficient  than  what  had  been  in  place. 

10.  The  routine  refrac  was  improved  to  proportionally  distribute  refracted  wave  energy  into 
adjacent  directional  bins.  Earlier  versions  of  this  routine  refracted  energy  into  one  bin  only  and 
thus  produced  and  unrealistic  step-like  patterns  in  directional  wave  spectra. 
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4.0  ARCHITECTURAL  DESIGN 

„CTI  V?  ^^ectund  Desi8“  section  shows  the  overall  design  and  the  calling  sequence  for  all 
Us  of  the  SURF.  Each  CSU  is  shown  in  the  calling  sequence  with  the  associated  CSU  related  to 
each  specific  unit.  Figure  (1)  presents  the  path  in  which  each  CSU  is  called  and  all  associated  CSUs 
which  in  turn  are  called  from  the  parent  unit.  Specific  details  concerning  the  criteria  for  each  CSU 
being  called  are  defined  in  the  Section  5.0:  SURF  CSU  Detailed  Design. 
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5.0 


CSCI  DETAILED  DESIGN 


5.1  Program  SURF 
Program  Call: 

SURFQ 

Summary: 


The  SURF  routine  is  the  starting  point  for  executing  SURF.  The  routine  identifies  the  input 
type  and  controls  the  reading  of  data  and  user  selected  computation  options.  The  routine  calls  the 
main  wave  parameter  calculation  routines  and  controls  the  output  of  the  resulting  data. 

Input  Variables:  None. 

Output  Variables:  None. 

Local  Variables: 


alfa 

Real 

bravo 

Real 

chrlie 

Real 

dangle 

Real 

depname 

Char*40 

dsea 

Real 

dstart 

Real 

dswell 

Real 

dxyl  (points) 

Real 

echo 

Real 

ehsig 

Real 

esowm  (dirNum.freqNum) 

Real 

file_dat 

Char*40 

file_in 

Char*40 

file_out 

Char*40 

file_tmp 

Char*40 

foxtrt 

Real 

fracname 

Char*40 

fireq  (freqNum) 

Real 

freql  (freqNum) 

Real 

freq2  (freqNum) 

Real 

Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Angle  Between  Directional  Bins 
Depth  Profile  File  Name 
Input  Direction  for  Sea  Contribution 
Input  Starting  Depth 

Input  Swell  Direction  for  Internally  Generated 
Spectrum 

Corresponding  Depths  with  No  Tide 
Breaker  Angle 

Significant  Wave  Height  from  Directional 
Spectrum 

Directional  Wave  Spectrum 
Output  File  Name  *.dat 
Input  Filename 
Output  File  Name  *.out 
Temporary  File 

Longshore  Current  Speed  and  Direction 
Wave  Refraction  File  Name 
Input  Wave  Spectrum  Center  Frequencies 
Beginning  Frequency  Bin  Values 
Ending  Frequency  Bin  Values 
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gamma2 

Real 

golfl 

Real 

golf? 

Real 

gt_frq 

Integer 

hsea 

Real 

hswell 

Real 

iday 

Integer 

idirec 

Integer 

iffeq 

Integer 

igamma 

Integer 

ihour 

Integer 

ihtll 

Real 

ihtl2 

Real 

imin 

Integer 

imonth 

Integer 

iyear 

Integer 

jgamma 

Integer 

line 

Char*80 

lin_stress 

Logical 

lndname 

Char*40 

nnn 

Integer 

pet  (4) 

Real 

period  (freqNum) 

Real 

psea 

Real 

pswell 

Real 

roller 

Logical 

selfjst 

Char*l 

slope 

Real 

spectra 

Logical 

spefile 

Char*40 

surfy 

Logical 

tide 

Real 

Beach  Orientation,  Compass  Heading  Directly 

Toward  Beach 

Number  of  Surf  Lines 

Surf  Zone  Width 

Spectrum  Type 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Pierson  Moskowitz  Spectrum 
Input  Significant  Wave  Height  for  Internally 
Generated  Spectrum 
Input  Day 

Number  of  Direction  Bins  in  the  Input 
Spectrum 

Number  of  Frequency  Bands  in  the  Input 
Spectrum 

Beach  Orientation  Rotated  90°  from  Original 
Heading  Toward  Beach 
Input  Hour 

Wind  Speed  Coded  Surf  Forecast  Value 

Wind  Direction 

Input  Minute 

Input  Month 

Input  Year 

Temporary  Value  Set  to  Beach  Orientation 

Temporary  Character  Variable 

Longshore  Current  Solution  (True  or  False) 

Input  Landing  Zone  Name 

Number  of  Points  in  the  Input  Depth  Array 

Percent  of  Different  Breaker  Types 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Period  Array  (1/Frequency) 

Input  Wave  Period  for  Sea  Contribution  to 

Pierson  Moskowitz  Spectrum 

Input  Swell  Period  for  Internally  Generated 

Spectrum 

Roller  Usage  (True  or  False) 

Self  Start  Flag  (Yes  or  No) 

Bottom  Slope 

Does  Input  Spectra  Exist?  (True  or  False) 
Selected  Wave  Spectrum  File  Name 
Significant  Wave  Heights  Greater  than 
0.5  ft?  (True  or  False) 

Input  Tide  Level 
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wdir  Real 

wspd  Real 

xcoeff  (dirNum,  freqNum)  Real 
xdelt  Real 

xdelt_gr  Real 

xfrom  (dirNum)  Real 

xtheta  (dirNum, freqNum)  Real 
xx  1  (points)  Real 

ydepth  Char*l 

ydetail  Char*l 

yrefrac  Char*l 

ystr  Char*l 


Subroutines  Called  from  SURF  (  ): 

ABORT 

CALCSURF 

DEPDRVR 

GENSPEC 

MDSRF1 

PRT_OUT3 

READRFRC 

READSPEC 

REFRAC 

SRFSETUP 

STRFRAC 

SUMMARY 

SURFCAST 


Input  Wind  Direction  Compass  Heading  Wind 
Blows  from 
Input  Wind  Speed 

Wave  Height  Refraction  Coefficients 
Surf  Zone  Output  Interval 
Self  Adjusting  Cross-Shore  Grid  Step 
Direction  Array,  Direction  Wave  Energy 
Comes  From 

Angle  Refraction  Coefficients 
Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

Input  Depth  Profile  Used?  (Yes  or  No) 
Detailed  Output?  (Yes  or  No) 

Is  Refraction  Considered  in  Analysis? 

(Yes  or  No) 

Is  Straight  Coast  Refraction  Used? 

(Yes  or  No) 
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Figure  2.  Program  SURF  Flowchart 


\ 


\ 


/Write Additonal  j 
Information  to  j 
OolpqtFfle  j 

\ 


\ 


Can 

SURFCAST 


5.2  Subroutine  ABORT 


Subroutine  Call: 
ABORT ( ) 


Summary: 

Subroutine  ABORT  acts  as  the  single  program  termination  routine.  The  subroutine  handles 
normal  program  execution  and  error  interrupt.  ABORT  is  called  to  stop  the  execution  of  the 
program.  If  an  error  interrupt  calls  ABORT,  the  error  message  is  generated  locally  in  the  calling 


Input  Variables:  None. 

Output  Variables:  None. 

Local  Variables:  None. 

Subroutines  Called  from  ABORT  (  ):  None. 

ABORT  (  )  Called  from  Subroutines: 

C_IN_DEP 

EQUILPRFMAIN_WAV 

MDSRF2 

NEW.BRK 

NONL1N2 

PRT_OUTl 

PRT_OUT2 

READRFRC 

READSPEC 

SRFSETUP 

SURF 
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Figure  3.  Subroutine  ABORT  Flowchart 
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53  Subroutine  BALANCEQ 


Subroutine  Call: 

BALANCEQ  (roller,  theta,  Cg,  rhs,  hrmsl,  dp,  mean_freq,  xk,  hrxns2,  convg) 

Summary: 

Subroutine  BALANCEQ  computes  new  wave  height  value  at  the  next  onshore  grid  cell  by 
performing  an  iterative  solution  to  the  energy  equations. 


Input  Variables: 

Cg 

Real 

Wave  Group  Velocity 

dp 

Real 

Offshore  Water  Depth 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

mean_frq 

Real 

Wave  Frequency 

rhs 

Real 

Right  Hand  Side  of  Energy  Balance 

Equation 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

theta 

Real 

Wave  Angle 

xk 

Real 

Wave  Number 

Output  Variables: 

convg 

Logical 

Convergence  Flag  (True  or  False) 

hrms2 

Real 

Significant  Wave  Height  at  next 

Onshore  Grid 

Local  Variables: 

avgh 

Real 

Average  Wave  Height 

check 

Real 

Convergence  Check 

done 

Logical 

Flag  indicating  End  of  Loop 

f3 

Real 

Function  which  Calculates  Total  Energy 

kount 

Integer 

Loop  Iteration  Counter 

lhs 

Real 

Left  Hand  Side  of  the  Energy  Equation 

limit 

Logical 

Flag  for  Comparison  of  the  Left  &  Right  Side 
of  the  Energy  Equation  (True  or  False) 

lowerh 

Real 

Lower  Limit  of  Wave  Height 

maxkount 

Integer 

Maximum  Number  of  Loop  Iterations  =1000 

oldavgh 

Real 

Previous  Average  Wave  Height  Value 

pet 

Real 

Convergence  Step  Value 

tol 

Real 

Convergence  Tolerance 

upperh 

Real 

Upper  Limit  of  the  Wave  Height 
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Subroutines  Called  from  BALANCEQ  (  ):None 
Functions  Called  from  BALANCEQ  (  ): 

F3 

BALANCEQ  (  )  Called  from  Subroutines: 

MAIN_WAV 

SLFJSTRT 


figure  4.  Subroutine  BALANCEQ  Flowchart 

f  START  1 


I  Call  F3 


IsLHS  < 
RHS? 


5.4  Subroutine  C_F1NE 
Subroutine  Call: 


C_FINE  (ndepth,  xxin,  zzin,  xdelt_gr,  nnn,  xxl,  dxyl) 
Summary: 


Subroutine  CFINE  linearly  interpolates  the  input  water  depths  and  offshore  distances  to  an 

evenly  spaced  grid.  The  internally  defined  grid  self-adjusts  to  maximize  spatial  resolution  without 
exceeding  the  array  dimensions. 


Input  Variables: 

ndepth 

Integer 

xdelt_gr 

Real 

xxin  (points) 

Real 

zzin  (points) 

Real 

Output  Variables: 

dxyl  (points) 

Real 

nnn 

Integer 

xxl  (points) 

Real 

Local  Variables: 

dxl 

Real 

dx2 

Real 

dxx 

Real 

dzz 

Real 

mm 

Integer 

mml 

Integer 

mmm 

Integer 

nn 

Integer 

xlast 

Real 

xtemp 

Real 

Subroutines  Called 

from  C_FINE  ( ); 

Number  of  Points  in  Input  Depth  Profile 
Self-Adjusting  Cross-Shore  Grid  Step 
Cross-Shore  Distances 
Corresponding  Depths 


Corresponding  Depths  without  Tide 
Number  of  Points  in  the  Input  Depth  Array 
Adjusted  Cross-Shore  Distances  from 
Depth  Profile 


Temporary  Variable  Used  in  Calculation  of 

Next  Grid  Point  Distance 

Temporary  Variable  Used  in  Calculation  of 

Next  Grid  Point  Distance 

Distance  Quotient 

Difference  Between  Depth  and 

Previous  Depth 

Counter  Variable 

Counter  Variable 

Counter  Variable 

Counter  Variable 

Last  Distance  Offshore  from  Input  Profile 
Temporary  Variable  for  Cross-Shore  Values 

None. 
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C_FINE  (  )  Called  from  Subroutines: 


C_REGR1D 

Figure  5.  Subroutine  C_FINE  Flowchart 


5.5  Subroutine  C_GAMMA 
Subroutine  Call: 

C_GAMMA  (gamma2,  igamma) 

Sununary: 

Subroutine  C_GAMMA  rotates  the  beach  orientation  data  read  from  the  input  file.  The  user 
defines  the  beach  orientation  as  the  compass  heading  of  a  boat  traveling  directly  toward  the  shore  on 
a  perpendicular  line  to  the  coast.  The  value  is  then  rotated  to  reflect  the  orientation  of  the  local 
coastline  with  respect  to  magnetic  north. 

Input  Variables: 


gamma2 

Real 

Beach  Orientation,  Heading  Directly 
Toward  Beach 

Output  Variables: 

igamma 

Integer 

Rotated  Beach  Orientation 

Local  Variables: 

gammatp 

mtemp 

Real 

Integer 

Temporary  Variable  Used  in  Calculation 
Temporary  Variable  in  Calculation 

Subroutines  Called  from  C_GAMMA  ( ): 

None. 

C_GAMMA  ( )  Called  from  Subroutines: 


i  DEPDRVR 

!■ 

| 

Figure  6.  Subroutine  C_GAMMA  Flowchart 

i 

f  START  ") 

f.  r-.  I 

Rotation  of  Gamma 
Relative  to  the  Beach 

i;i  _ 

t  i 

[  RETURN  ^ 
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5.6  Subroutine  C_IN_DEP 
Subroutine  Call: 

C_IN_DEP  (depname,  spedepth,  xdelt_gr,  nnn,  xxl,  dxyl,  dstart) 

Summary: 

Subroutine  C_IN_DEP  reads  the  depth  profile  and  header  information  contained  in  the  input 
data  file.  The  routine  identifies  the  units  of  measurement  used  to  construct  the  depth  profile  and 
checks  the  order  of  the  offshore  distances.  If  the  data  is  misaligned,  the  subroutine  will  sort  and 
reorder  the  depths  and  offshore  distances. 


Input  Variables: 

depname 

Char*40 

Depth  Profile  File  Name 

dstart 

Real 

Input  Starting  Depth 

xdelt_gr 

Real 

Self  Adjusting  Cross-Shore  Grid  Step 

Output  Variables: 

dxyl  (points) 

Real 

Corresponding  Depths  without  Tide 

nnn 

Integer 

Number  of  Points  in  the  Input  Depth  Array 

xxl  (points) 

Real 

Adjusted  Cross-Shore  Distances  from  the 
Depth  Profile 

Local  Variables: 

al 

Real 

Temporary  Variable 

a2 

Real 

Temporary  Variable 

adum 

Char*80 

Temporary  Variable,  Character  String  in  Input 
Field 

dcall 

Real 

Conversion  Factor  for  Distance  Offshore, 
Convert  to  Meters 

dcal2 

Real 

Conversion  Factor  for  Depths  Offshore, 
Convert  to  Meters 

dx 

Real 

Temporary  Variable  for  Distance  Offshore 
from  Input  File 

dz 

Real 

Temporary  Variable  for  Depths 

I 

Integer 

Loop  Variables 

icall 

Integer 

Input  from  Depth  File, 

Units  of  Distance  Offshore 

1  =  Feet 

*cal2  Integer 


instat 

Integer 

j 

Integer 

k 

Integer 

line 

Integer 

loop 

Integer 

ndepth 

Integer 

xxin  (points) 

Real 

zzin  (points) 

Real 

Subroutines  Called  from  C_IN_DEP  ( ): 

ABORT 

C_UN 

C_REGRID 

C_IN_DEP  (  )  Called  from  Subroutines: 
DEPDRVR 


2  =  Meters 

3  =  Yards 

Depth  Units  Input  from  Depth  File 

1  =  Feet 

2  =  Meters 

3  =  Fathoms 
Hie  Open  Status 
Loop  Variables 

Temporary  Variable  for  Number  of  Points 
Counter  for  the  Number  of  Lines  in  the  Input 
Depth  Profile 
Loop  Counter 

Number  of  Points  in  Input  Depth  Profile 
Cross-Shore  Distances 
Corresponding  Depths 
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Figure  7.  Subroutine  CJNJDEP  Flowchart 


[  START  } 


\ 


Set  Unit 
Conversion 
Factors 


J 


Decrease 
Starting  Depth 


t 


End  of  Loop 


I 

/ Find  Array  Index  ^ 
Where  Depth  is 
Approximately  Equal 
to  the  Starting  Depth 


Call 

C_REGRID 


c 


RETURN 


j 


Increase 
Starting  Depth 


Call  ABORT 


Print  Warning 
Message 
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5.7  Subroutine  C.REGRID 
Subroutine  Call: 


C.REGRID  (ndepth,  xxin,  zzin,  xdelt_.gr,  nnn,  xxl,  dxyl) 

Summary: 

Subroutine  C.REGRID  examines  the  cross-shore  step  size  (Ax)  of  the  input  depth  profile  and 
selects  a  new  step  size  to  optimize  the  depth  and  cross-shore  distance  arrays.  The  step  size  is 
automatically  adjusted  and  the  arrays  are  constructed  so  the  length  does  not  exceed  the  dimension 
of  the  array. 


Input  Variables: 


ndepth  Integer 

xdelt_gr  Real 

xxin  (points)  Real 

zzin  (points)  Real 

Output  Variables: 

nnn  Integer 

xdelt.gr  Real 

xxl  (points)  Real 

xxin  (points)  Real 

zzin  (points)  Real 

Local  Variables:  None. 


Subroutines  Called  from  C.REGRID  ( ): 
C.FINE 

C.REGRID  ( )  Called  from  Subroutines: 
C.IN.DEP 


Number  of  Points  in  Depth  Profile 
Self  Adjusting  Cross-Shore  Grid  Step 
Cross-Shore  Distances 
Corresponding  Depths 


Number  of  Points  in  Input  Depth  Array 
Self  Adjusting  Cross-Shore  Grid  Step 
Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

Adjusted  Cross-Shore  Distances 
Corresponding  Depths 
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Figure  8.  Subroutine  C_REGRID  Flowchart 


[  START  ] 

1 

S  Number  of  ^ 
Points  in  Input 
Depth  Array 
_ Loop _ 


\ 


f  RETURN  j 
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5.8  Subroutine  CUN 


Subroutine  Call: 

C_UN  (dcall,  dcal2,  ndepth,  xxin,  zzin,  xdelt_gr,  spedepth) 

Summary: 

Subroutine  C_UN  converts  measurement  units  of  input  cross-shore  distances,  depth  arrays, 
starting  depth  and  the  grid  step  size  (Ax)  to  meters  for  internal  calculations. 


Input  Variables: 

dcall 

Real 

dcal2 

Real 

dstart 

Real 

ndepth 

Integer 

xdelt_gr 

Real 

xxin  (points) 

Real 

zzin  (points) 

Real 

Output  Variables: 

spedepth 

Real 

xdelt _gr 

Real 

xxin  (points) 

Real 

zzin  (points) 

Real 

Local  Variables: 

I  Integer 

Subroutines  Called  from  C_UN  ( ): 
C_UN  ( )  Called  from  Subroutines: 


Conversion  Factor  for  Cross-Shore  Distances 

Conversion  Factor  for  Water  Depths 

Input  Starting  Depth 

Number  of  Points  in  Input  Depth  Profile 

Self  Adjusting  Cross-Shore  Grid  Step 

Cross-Shore  Distances 

Corresponding  Depths 


Input  Starting  Depth 
Self  Adjusting  Cross-Shore  Grid  Step 
Cross-Shore  Distances 
Corresponding  Depths 


Loop  Counter 
None. 


Figure  9.  Subroutine  C_UN  Flowchart 
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5.9  Subroutine  CALC_HB3 


Subroutine  Call: 

CALC_HB3  (dp,  hrms,  p_flag,  hb3) 

Summary: 

Subroutine  CALC_HB3  integrates  the  wave  height  distribution  for  a  given  toot  mean  square 
wave  height  and  calculates  a  term  inherent  to  the  roller  dissipation  function. 

Input  Variables: 

dP  Real  Offshore  Water  Depth 

hr®s  Reaj  Root  Mean  Square  Wave  Height  Calculation 

P-S  Logical  Weighting  Factor  Flag  (True  or  False) 

Output  Variables: 


hb3  Real 

Local  Variables: 


Weighting  Function  for  Dissipation  Term 


Real  Maximum  Wave  Height 

[dow  Real  Minimum  Wave  Height 

integrat  Real  Wave  Height  Distribution  Calculated  for  a 

Single  Wave  at  a  Specific  Location 

Functions  Called  from  CALC_HB3  ( ): 

INTEGRAT 

CALCJHB3  (  )  Called  from  Subroutines: 


CALCROLL 

GET_DISS 
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Figure  10.  Subroutine  CALC_HB3  Flowchart 


31 


5.10  Subroutine  CALCROLL 
Subroutine  Call: 

CALCROLL  (roller,  hrms,  dp,  fqz,  theta,  xk,  e_roller) 

Sununary: 

Subroutine  CALCROLL  calculates  roller  energy  at  a  point  in  the  surf  zone  based  on  water 
depth  and  Wave  Height  (hrms)  at  that  location. 


Input  Variables: 

dp 

Real 

Offshore  Water  Depth 

rqz 

hrms 

roller 

theta 

Real 

Zero  Crossing  Frequency 

Real 

Root  Mean  Square  Wave  Height 

Logical 

Roller  Option  Flag  (True  or  False) 

Real 

Wave  Angle,  Representative  of  Radiation 

xk 

Real 

Stress  Angle 

Wave  Number 

Output  Variables: 

e_roller 

Real 

Roller  Contribution  to  Energy  Equation 

Local  Variables: 

c 

Real 

Wave  Celerity 

er 

hb3 

Real 

Temporary  Roller  Variable 

Real 

Weighting  Function  for  Dissipation  Term 

P_flag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

z 

Real 

Roller  Energy  Multiplier 

Subroutines  Called 

from  CALCROLL( ): 

CALC_HB3 


CALCROLL(  )  Called  from  Subroutines: 
GET_RHS 
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CALCROLL(  )  Called  from  Functions: 


Figure  12.  Subroutine  CALCROLL  Flowchart 
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5.11  Subroutine  CALCSURF 


Subroutine  Call: 

,  CALCSURF  ( roller,  lin.stress,  ehsig,  wspd,  wdir,  tide,  ydepth,  nnn,  dxyl,  xxl,  ifieq,  freql, 

neq2,  fteq,  idirec,  xfirom,  esowm,  dstart,  igamma,  ydetail,  iyear,  imonth,  iday,  ihour,  imin,  xdelt, 

xdelt_gr,  self_st,  file_spc,  surf,  pet,  alfa,  bravo,  chrlie,  echo,  foxtrt,  golfl,  golf2,  ihtll  ihtl2 
jgamma) 

Sununary: 

Subroutine  CALCSURF  acts  as  the  primary  driver  for  the  various  subroutines,  which 
calculate  wave  parameters  and  the  longshore  current  across  the  surf  zone. 


Input  Variables: 

dstart 

Real 

dxyl  (points) 

Real 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

freq  (freqNum) 

Real 

freql  (freqNum) 

Real 

freq2  (freqNum) 

Real 

iday 

Integer 

idirec 

Integer 

ifreq 

Integer 

igamma 

Integer 

ihour 

Integer 

imin 

Integer 

imonth 

Integer 

iyear 

Integer 

lin_stress 

Logical 

nnn 

Integer 

roller 

Logical 

selfjst 

Char*l 

tide 

Real 

wdir 

Real 

wspd 

Real 

xdelt 

Real 

xdelt_gr 

Real 

xfrom  (dirNum) 

Real 

Input  Starting  Depth 
Corresponding  Depths  without  Tide 
Significant  Wave  Height  from 
Directional  Spectrum 
Directional  Wave  Spectrum 
Input  Wave  Spectrum  Center  Frequency 
Beginning  Frequency  Bin  Value 
Ending  Frequency  Bin  Value 
Input  Day 

Number  of  Direction  Bins  in  Input  Spectrum 
Number  of  Frequencies  in  Input  Spectrum 
Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 
Input  Hour 
Input  Minute 
Input  Month 
Input  Year 

Longshore  Current  Solution  (True  or  False) 
Number  of  Points  in  Input  Depth  Array 
Roller  Option  Flag  (True  or  False) 

Self  Start  Flag  (Yes  or  No) 

Input  Tide  Level 

Input  Wind  Direction,  Compass  Heading 

Wind  is  Blowing  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Self-Adjusting  Cross-Shore  Grid  Step 

Direction  Array,  Direction  Wave  Energy 
Comes  From 
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xx  1  (points) 

Real 

ydepth 

Char*  1 

ydetail 

Char*l 

Output  Variables: 

alfa 

Real 

bravo 

Real 

chrlie 

Real 

echo 

Real 

foxtrt 

Real 

golfl 

Real 

golf2 

Real 

ihtll 

Real 

ihtl2 

Real 

jgamma 

Integer 

pet  (4) 

Real 

surf 

Logical 

Local  Variables: 

along  (points) 

Real 

b 

Real 

bl  (points) 

Real 

blong  (points) 

Real 

c 

Real 

cl 

Real 

c2 

Real 

c3 

Real 

c4 

Real 

cf 

Real 

Cg 

Real 

clong  (points) 

Real 

convg 

Logical 

df 

Real 

Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

Input  Depth  Profile  Used?  (Yes  or  No) 
Detailed  Output?  (Yes  or  No) 


Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Breaker  Angle 

Longshore  Current  Speed  and  Direction 

Number  of  Surf  Lines 

Surf  Zone  Width 

Wind  Speed 

Wind  Direction 

Temporary  Value  Set  to  Beach  Orientation 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Flag  for  Low/No  Surf  Conditions 
(True  or  False) 


Horizontal  Mixing  Parameter  from 
Thornton  &  Whittord 

Empirical  Factor  in  Thornton  &  Guza  Wave 
Breaking  Model  (=  1.00) 

Bottom  Slope 

Bottom  Friction  for  Deep  &  Shallow  Water 

Wave  Celerity  at  Input  Starting  Depth 

Eddy  Viscosity  Coefficient 

Bottom  Friction  Coefficient 

Radiation  Stress  Coefficient  □  Multiple  for 

Longshore  Current  Model 

Longshore  Wind  Stress  Coefficient  □  Multiple 

for  Longshore  Current  Model 

Coefficient  of  Bottom  Friction 

Wave  Group  Velocity 

Wind  Stress  Contribution  to 

Longshore  Current 

Energy  Equation  Convergence  Flag 

Difference  Between  Adjacent  Frequency  Bins 
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distmax 

Real 

dp 

Real 

dth 

Real 

dws_stop 

Integer 

dxy  (points) 

Real 

eb_last 

Real 

ebtemp  (points) 

Real 

file_spc 

Char*40 

fqd 

Real 

fqz 

Real 

ftsq2msq 

Real 

hlmax 

Real 

h2max 

Real 

hrms 

Real 

htemp  (points) 

Real 

iimax 

Integer 

irealf 

Integer 

j 

Real 

jJi 

Integer 

j_ii2 

Integer 

k 

Real 

per 

Real 

print_spc 

Integer 

ptemp  (points) 

Real 

ik  (points,4) 

Real 

stringout 

Character 

stringsub 

Character 

suml 

Real 

temp 

Real 

theta 

Real 

thetal 

Real 

Farthest  Offshore  Distance 
Offshore  Water  Depth 

Difference  Between  Adjacent  Directional  Bins 
Flag  for  Shallow  Water  Directional  Wave 
Spectrum  Print  Control 
Pre-Tidal  Depth  with  Tide 
Roller  Dissipation  Term  at  Farthest 
Point  Offshore 

Temporary  Roller  Dissipation  Term 
Across  Transect 

File  Name  of  Shallow  Water  Directional  Wave 
Spectrum 

Peak  Frequency  at  the  Center  of  the  Frequency 
Band 

Zero  Crossing  Frequency 

Conversion  Factor  from  Feet  Squared  to  Meters 

Squared 

Largest  Significant  Wave  Height  in  the 
Surf  Zone 

Largest  Maximum  Wave  Height  in  the 
Surf  Zone 

Root  Mean  Square  Wave  Height 
Temporary  Variable  for  Significant  Wave 
Height  Values 

Number  of  Calculation  Locations 
Cutoff  Index  for  Printing  Shallow  Water 
Directional  Wave  Spectrum 
Temporary  Variable  for  Cross-Shore  Values 
Index  where  Wave  Probabilities  come 
Above  Threshold 

Longshore  Current  Loop  Variable  for  Outer 
Edge  of  Surf  Zone 
Temporary  Variable  for  Significant 
Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 
Flag  to  Print  Shallow  Water  Wave  Spectrum 
Percentage  of  Breaking  Waves  and 
Breaker  Types 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 

Shallow  Water  Wave  Spectrum  Output  String 

Temporary  String  Variable 

Sum  of  Wave  Length  in  the  Surf  Zone 

Temporary  Variable 

Wave  Angle 

Wave  Angle  at  Input  Starting  Depth 
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theta2  Real 

v  (points)  Real 

vraax  Real 

vmin  Real 

vwind  Real 

wdspd  Real 

wid_ii  Integer 

width  Real 

xk  Real 

xktemp  (points)  Real 

xshift  Real 

xtemp  (points)  Real 


Subroutines  Called  from  CALCSURF  (  ): 

GT_SIG_H 

INTTLIZE 

KLONG 

MAIN_WAV 

NONLIN 

PRT_OUTl 

PRT_OUT2 

RAD_ST1 

S_COEFF 

S_NOSURF 

S_TIDE 

SHORTOUT 

ZONE1 

CALCSURF  (  )  Called  from  Subroutines: 
SURF 


Wave  Angle  at  Input  Starting  Depth 
Longshore  Current 

Maximum  Positive  Longshore  Current 

Maximum  Negative  Longshore  Current 

Group  Wind  Velocity 

Wind  Speed  Conversion 

Knots  to  CM/S  =  51.44 

Array  Location  for  Surf  Zone  Width 

Surf  Zone  Width 

Wave  Number 

Temporary  Variable  for  Wave  Number 
Horizontal  Cross-Shore  Location 
Temporary  Variable  for  Cross-Shore  Values 
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5.12  Subroutine  DEPDRVR 


Subroutine  Call: 

DEPDRVR  (depname,  spedepth,  xdelt,  ydepth,  slope,  gamma2,  nnn,  xxl,  dxyl,  igamma, 
xdelt_gr,  dstart,  ystr) 

Summary: 

Subroutine  DEPDRVR  is  the  driver  routine  for  reconstructing  the  depth  arrays  in  an  optimized 
step  size. 


Input  Variables: 


depname  Char*40 

dstart-  Real 

gamma2  Real 

slope  Real 

xdelt  Real 

ydepth  Char*l 

ystr  Char*l 

Output  Variables: 

dxyl  (points)  Real 

igamma  Integer 

nnn  Integer 

xdelt_gr  Real 

xxl  (points)  Real 

Local  Variables:  None. 


Subroutines  Called  from  DEPDRVR  (  ): 

EQUILPRF 

CjGAMMA 

CJN.DEP 


DEPDRVR  ( )Called  from  Subroutines: 
SURF 


Depth  Profile  File  Name 

Input  Starting  Depth 

Beach  Orientation  Compass  Heading 

Directly  Toward  Beach 

Bottom  Slope 

Surf  Zone  Output  Interval 

Usage  of  Input  Depth  Profile  (Yes  or  No) 

Straight  coast  refraction  flag 


Corresponding  Depths  without  Tide 
Beach  Orientation  Rotated  90  Degrees  from  the 
Original  Heading  Toward  the  Beach 
Number  of  Points  in  the  Input  Depth  Array 
Self-Adjusting  Cross-Shore  Grid  Step 
Adjusted  Cross-Shore  Distances  from  the 
Depth  Profile 
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Figure  13.  Subroutine  DEPDRVR  Flowchart 


5.13  Subroutine  EQUILPRF 


Subroutine  Call: 

EQUILPRF  (wavdep,  ystr,  rtype,  xgrd,  numstep,  xxl,  dxyl,  deepest_depth) 

Sununary: 


Subroutine  EQUILPRF  constructs  a  depth  profile  for  surf  calculations.  This  equilibrium 
profile  is  based  on  the  equation:  y=AxA(2/3),  where  A  is  a  coefficient  related  to  sediment  grain  size  or 
frictional  dissipation.  This  equation  was  developed  by  Dean  (1977)  from  a  study  of  more  than  200 
beach  profiles.  The  "A"  coefficient  in  the  equilibrium  equation  has  units  of  meters,  calculations  in 
feet  require  different  values  or  conversion  to  feet  after  initial  calculations.  Sediment/grain  types  are 
denoted  by  the  variable  "rtype"  which  is  the  index  for  a  value  in  the  array  of  coefficients  defining  the 
following  grain  sizes: 

1  =  boulders 

2  =  cobble 

3  =  pebbles 

4  =  granules 

5  =  very  coarse  sand 

6  =  coarse  sand 

7  =  medium  sand 

8  =  fine  sand 

9  =  very  fine  sand 

10  =  silt 


Input  Variables: 

dpthoff 

Real 

Input  Starting  Depth 

numstep 

Integer 

Number  of  Points  in  the  Input  Depth  Array 

rtype 

Real 

Sediment/  Grain  Type 

wavdep 

Real 

Input  Wave  Depth 

xgrd 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

ystr 

Character 

Straight  Coast  Refraction  Logic 

Output  Variables: 

dxyl  (points) 

Real 

Corresponding  Depths  with  No  Tide 

xxl  (points) 

Local  Variables: 

Real 

Cross-Shore  Distances 

a(10) 

Real 

Array  of  Sediment  Coefficients 

ause 

Real 

Actual  Sediment  Type  Coefficient  for  Profile 

call 

Real 

Conversion  Factor  (Meters) 
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distance 

diston 

dpthon 

I 

x 

xone 

z 


Logical 

Real 

Real 

Integer 

Real 

Real 

Real 


Subroutines  Called  from  EQUILPRF  ( ): 
EQUILPRF  ( )  Called  from  Subroutines: 
DEPDRVR 


Flag  for  Equilibrium  Depth  Bottom 
Highest  Onshore  Distance 
Highest  Onshore  Depth 
Loop  Counter 
Temporary  Variable 
Farthest  Point  Offshore 
Temporary  Variable 

ABORT 
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Figure  14.  Subroutine  EQUILPRF  Flowchart 


5.14  Subroutine  GENSPEC 


Subroutine  Call: 

GENSPEC  (hsea,  psea,  dsea,  hswell,  pswell,  dswell,  ifreq,  idirec,  freq,  freql,  freq2,  xfrom, 
esowm,  period,  ehsig,  dangle) 

Sununary: 

Subroutine  GENSPEC  initializes  matrices  for  the  creation  of  an  internally  generated 
directional  wave  spectrum.  This  wave  spectrum  has  50  frequencies  and  36  directions. 

Input  Variables: 


dsea 

Real 

dswell 

Real 

hsea 

Real 

hswell 

Real 

psea 

Real 

pswell 

Real 

Output  Variables: 

dangle 

Real 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

freq  (freqNum) 

Real 

freql  (freqNum) 

Real 

freq2  (freqNum) 

Real 

idirec 

Integer 

ifreq 

Integer 

period  (freqNum) 

Real 

xfrom  (dirNum) 

Real 

Input  Direction  for  Sea  Contribution 

Input  Swell  Direction  for  Internally 

Generated  Spectrum 

Input  Significant  Wave  Height  for  Sea 

Contribution  to  Pierson  Moskowitz  Equation 

Input  Significant  Wave  Height  for  Internally 

Generated  Spectrum 

Input  Wave  Period  for  Sea 

Contribution  to  Pierson  Moskowitz  Equation 

Input  Swell  Period  for  Internally 

Generated  Spectrum 


Angle  Between  Directional  Bins 

Significant  Wave  Height  from 

Directional  Spectrum 

Directional  Spectrum 

Input  Wave  Spectrum  Center  Frequencies 

Beginning  Frequency  Bin  Values 

Ending  Frequency  Bin  Values 

Number  of  Direction  Bins  in  the 

Input  Spectrum 

Number  of  Frequencies  in  the 

Input  Spectrum 

Period  Array  (1/Frequency) 

Direction  Array,  Direction  Wave  Energy 
Comes  From 
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Local  Variables: 

df  Real  Difference  between  Frequency  Bins 

idir  Integer  Direction  Loop  Counter 

ifirq  Integer  Frequency  Loop  Counter 

Subroutines  Called  from  GENSPEC  ( ): 

WAVEFTT  * 

GENSPEC  (  )  Called  from  Subroutines: 

SURF 

Figure  15.  Subroutine  GENSPEC  Flowchart 
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5.15  Subroutine  GET_BRK 


Subroutine  Call: 

GET_BRK  ( ii,  nnn,  xxl,  dxy,  xdelt_gr,  hrms,  10,  per,  xoff,  rk,  bl,  brklO,  distmax,  p) 
Summary: 

Subroutine  GET_BRK  calculates  percentage  of  breakers  and  percent  breaker  type  given  at 
each  point  along  the  transect:  p  (1)  =  Spilling,  p  (2)  =  Plunging,  p  (3)  =  Surging,  p  (4)  =  100*Sum. 


Input  Variables: 

bl  (points) 

Real 

brklO 

Logical 

distmax 

Real 

dxy  (points) 

Real 

hrms 

Real 

ii 

Integer 

per 

Real 

rk  (points,4) 

Real 

xdelt_gr 

Real 

xoff 

Real 

Output  Variables: 

bl  (points) 

Real 

brklO 

Logical 

distmax 

Real 

P(4) 

Real 

rk  (points,4) 

Real 

Local  Variables: 

beta 

Real 

Bottom  Slope 

Flag  for  First  Location  where  10%  of  the 
Waves  are  Breaking  (True  or  False) 

Farthest  Offshore  Distance 

Adjusted  Depths  with  Tide 

Root  Mean  Square  Wave  Height 

Index  where  Wave  Probabilities  Exceed 

Threshold 

Peak  Period  of  Directional  Wave  Spectrum 
Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 
Self-Adjusting  Cross-Shore  Step 
Distance  Offshore 


Bottom  Slope 

Rag  for  First  Location  where  10%  of  the 
Waves  are  Breaking  (True  or  False) 

Farthest  Offshore  Distance 
Temporary  Array  for  Breaker 
Percentage  Totals 
Percent  Breaker  of  Each  Type 

Temporary  Variable  for  Bottom  Slope 
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Subroutines  Called  from  GET_BRK  (  ): 
PERCENT 

GET_BRK  (  )  CaUed  from  Subroutines: 
MAIN.WAV 
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5.16  Subroutine  GET_DISS 


Subroutine  Call: 

GET_DISS  (roller,  fqz,  dp,  hrms,  p_flag,  diss) 

Summary: 

Subroutine  GET_DISS  returns  the  wave  dissipation  factor.  This  term  is  based  on  a  bore 
dissipation  model  and  can  include  roller  dissipation  if  selected.  The  dissipation  term  is  included  in  the 
pwave  energy  balance  equation.  The  wave  dissipation  is  given  by: 

£  IPJJJe.  H3nm*M  * B3 


where  p  is  density,  g  is  gravity,  f  is  bottom  friction,  h  is  the  water  depth,  M  is  a  weighting  function 
based  on  hrms,  and  B  is  an  empirical  factor. 


Input  Variables: 


b 

Real 

Empirical  Factor  in  Thornton  &  Guza  Wave 
Breaking  Model  =1.00 

dp 

Real 

Offshore  Water  Depth 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

P_flag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

Output  Variables: 

• 

diss 

Real 

Bore  or  Roller  Dissipation  Function 
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Local  Variables: 


Real 

Real 


Weighting  Function  for  Dissipation  Term 
Dissipation  Function 


hb3 
z 

Subroutines  Called  from  GET_DISS  (  ): 
CALCHB3 

GET_DISS  (  )  Called  from  Subroutines: 
GET.RHS 

Figure  17.  Subroutine  GET_DISS  Flowchart 


I 


RETURN 


5.17  Subroutine  GET_M 


Subroutine  Call: 

GET_M  (dp,  hrms,  m) 

Summary: 

Subroutine  GET_M  calculates  the  weighting  function  multiplier. 

Input  Variables: 

dp  Real  Offshore  Water  Depth 

hrms  Real  Root  Mean  Square  Wave  Height 

Output  Variables: 

m  Real  Multiplier 

Local  Variables:  None. 

Subroutines  Called  from  GET_M  ( ):  None. 

GETJVf  (  )  Called  from  Subroutines: 

WEIGHTFN 

Figure  18.  Subroutine  GET_M  Flowchart 

f  START  ] 

I 

Calculate  Multiplier 
Based  Upon 
Weighting  Function 

i 

[  RETURN  ] 
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5.18  Subroutine  GETJP 


Subroutine  Call: 

GET_P  (frac,  p) 

Summary: 

Subroutine  GET_P  calculates  the  percentage  of  each  breaker  type  and  fills  the  conesponding 
array  elements. 

Input  Variables: 


frac  (3)  Real 

Output  Variables: 

P  (4)  Real 

Local  Variables: 

sum  Real 

Subroutines  Called  from  GET_P  ( ): 
GET_P  (  )  Called  from  Subroutines: 
PERCENT 


Temporary  Array  for  Breaker 
Percentage  Totals 


Percent  of  Different  Breaker  Types 
p  (1)  =  Spilling 
p  (2)  =  Plunging 
p  (3)  =  Surging 
p  (4)  =  Total 


Temporary  Variable  for  Total  of 
Percentage  Breakers 

None. 
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Figure  19.  Subroutine  GET_P  Flowchart 


[  START  ) 


1 

f  RETURN  ) 
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5.19  Subroutine  GET_RHS 
Subroutine  Call: 


GET_RHS  (roller,  hrms,  theta,  Cg,  dp,  xk,  b,  fqz,  xdelt_gr,  rhs,  diss) 


Summary: 


Subroutine  GET_RHS  calculates  the  right  hand  side  of  the  wave  energy  equation. 


Input  Variables: 
b 

Cg 

dp 

fqz 

hrms 

roller 

theta 

xdelt _gr 
xk 

Output  Variables: 

diss 

rhs 

Local  Variables: 

e_roller 

e_wave 

P_flag 


Real 

Real 

Real 

Real 

Real 

Logical 

Real 


Real 

Real 


Real 

Real 


Real 

Real 

Logical 


Subroutines  Called  from  GET_RHS  ( ): 


Empirical  Factor  in  Breaking  Model  =1.0 

Wave  Group  Velocity 

Offshore  Water  Depth 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Roller  Option  Flag  (True  or  False) 

Wave  Angle,  Representative  of  Radiation 
Stress  Angle 

Self-Adjusting  Cross-Shore  Grid  Step 
Wave  Number 


Bore  or  Roller  Dissipation  Function 
Right  Hand  Side  of  Energy  Equation 


Roller  Contribution  to  the  Energy  Equation 
Wave  Contribution  to  the  Energy  Equation 
Weighting  Factor  Flag  (True  or  False) 


CALCROLL 

GET_DISS 

GETJWAVE 


GET_RHS  ( )  Called  from  Subroutines: 


MAINJVAV 
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Figure  20.  Subroutine  GET_RHS  Flowchart 
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5.20  Subroutine  GETJSLOPE 
Subroutine  Call: 

GET_SLOPE  (h,  x,  iO,  lamda,  nnn,  beta) 

Summary: 

Subroutine  GET_SLOPE  gets  bottom  slope  from  a  depth  profile  for  percent  breaker 
calculations. 


Input  Variables: 


h 

Real 

Array  of  depths  associated  with  x  (meters) 

lamda 

Real 

Wavelength  computed  from  wavenumO 

X 

Real 

Array  of  cross-shore  distance  (meters) 

iO 

Integer 

index  from  x  where  slope  is  computed,  x(i0) 

nnn 

Integer 

number  of  elements  of  arrays  h  and  x  with  real  I 
nfo 

Output  Variables: 

Beta 

real 

bottom  slope;  positive  up  looking  at  beach 

Local  Variables: 


I  integer 

I_last_wet  integer 

Linshore  integer 

I_offshore  integer 

Once_dry  logical 

Slope_length  real 

Sum_of_slopes  real 

Temp  real 


Subroutines  Called  from  GET_SLOPE  (  ): 
none 


counter  in  do  loops 

index  of  last  underwater  in  h(x)  array 

index  of  x  and  h  used  for  slope  computation 

index  of  x  and  h  used  for  slope  computation 

used  in  getting  i_last_wet 

lamda/2.0 

temp  sum  of  slopes 

temporary  variable 


GET_RHS  ( )  Called  from  Subroutines: 

SLFJSTRT 

GET.BRK 
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5.21  Subroutine  GET_WAVE 
Subroutine  Call: 

GET_WAVE  (hrms,  theta,  Cg,  e_wave) 

Summary: 

Subroutine  GET_WAVE  calculates  wave  energy  flux  using  linear  wave  theory.  The  wave 
energy  flux  is: 


e=pamL 

8 


Cg  cos  0 


where  p  is  water  density,  g  is  gravity,  H  is  wave  height,  Cg  if  group  velocity,  and  0  is  the  wave 
angle. 


Input  Variables: 

Cg  Real 

hrms  Real 

theta  Real 

Output  Variables: 

e_wave  Real 

Local  Variables: 

ew  Real 

Subroutines  Called  from  GET_WAVE  (  ): 

GET_WAVE  (  )  Called  from  Subroutines: 


Wave  Group  Velocity 

Root  Mean  Square  Wave  Height 

Wave  Angle 


Energy  Flux 


Wave  Energy 
None. 


Figure  22.  Subroutine  GET_WAVE  Flowchart 
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5.22  Subroutine  GRID_FRC 
Subroutine  Call: 


GRID_FRC  (s,  r,  ifreq,  f,  idirec,  d,  ifreq_i,  fi,  idirecj,  di,  si,  ri) 

Summary: 

Subroutine  GRID_FRC  uses  a  bilinear  interpolation  scheme  to  re-grid  the  refraction  and 
shoaling  matrices  obtained  from  subroutine  readfrc  to  the  same  frequencies  and  directions  associated 
with  the  directional  wave  spectrum  obtained  from  the  subroutine  readspec. 


Input  Variables: 


d 

di 

Real 

Real 

direction  array  associated  with  input  refraction  and 
shoaling  matrices  (deg  +CW  from  N) 

direction  array  associated  with  wave  spectrum,  array  to 
which  input  refraction  and  shoaling  matrices  are 

interpolated  (deg  +CW  from  N) 

f 

fi 

Real 

frequency  array  associated  with  input  refraction  and 
shoaling  matrices  (hertz) 

Real 

integer 

frequency  array  to  which  input  refraction  and  shoaling 
matrices  are  interpolated  (hertz) 

idirec 

no.  of  direction  bands  in  the  input 
refraction  and  shoaling  matrices 

idirecj 

integer 

no.  of  direction  bands  to  which  the 
refraction  and  shoaling  fields  will  be 
interpolated 

ifreq 

integer 

no.  of  frequency  bands  in  the  original 
refraction  and  shoaling  matrices 

ifreqj 

integer 

no.  of  frequency  bands  to  which  the 
input  refraction  and  shoaling  matrices  will 
be  interpolated 

r(dimum,  freqnum) 

real 

input  refraction  matrix,  (units  are  degrees 
+CW  from  N.,  energy  from) 

s(dimum,  freqnum) 

real 

input  shoaling  matrix  (unitless;  energy 
shoaling  [kA2]) 

Output  variables: 

ri(dimum,  freqnum) 

real 

refraction  matrix  with  same  freqs  and  directions 
as  input  directional  wave  spectrum.(units  are  degrees 
(+CW  from  N.,  energy  from) 

si(dimum,  freqnum) 

real 

shoaling  matrix  with  the  same  freqs  and  directions 
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Local  Variables: 


i 

integer 

j 

integer 

k 

integer 

m 

integer 

tl 

real 

t2 

real 

Subroutines  Called  from  GRIDFRC  ( ) 
none 

GRID_FRC  ( )  Called  from  Subroutines 

Main  program  surf 


as  tiie  input  directional  wave  spectrum, 
(unitless;  energy  shoaling  [kA2]) 


counter 

counter 

counter 

counter 

temp 

temp 


61 


Figure  23.  Subroutine  GRIDJFRC  Flowchart 
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5.23  Subroutine  GRIDOUT 


Subroutine  Call: 

GRIDOUT  (ii,  xoffl,  xtemp,  dxy,  htemp,  ptemp,  xktemp,  thetatemp,  v,  dpi,  houtl,  hmax, 
pbreak,  brkrang,  vlngl) 

Summary: 

Subroutine  GRIDOUT  linearly  interpolates  parameters  for  final  output  using  the  user  defined 


cross-shore  step  width. 

Input  Variables: 

dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 

ii 

Integer 

Height  Values 

Index  where  Wave  Probabilities  Exceed 

ptemp  (points) 

Real 

Threshold 

Percentage  of  Breaking  Waves  and 

xktemp 

Real 

Breaker  Types 

Temporary  Variable  for  Wave  Number 

xoffl 

Real 

Distance  Offshore 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

v  (points) 

Real 

Longshore  Current 

Output  Variables: 

dpi 

Real 

Offshore  Depth 

hmax 

Real 

Maximum  Wave  Height  /  10.0 

houtl 

Real 

Significant  Wave  Height 

pbreak 

Real 

Percentage  Breaking  Waves 

vlngl 

Real 

Longshore  Current  Velocity 

wlen 

Real 

Wave  Length 

Local  Variables: 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

Subroutines  Called  from  GRIDOUT  ( ): 


GT_P 

UN_1 

UN_2 

LIN_3 

GRIDOUT  (  )  Called  from  Subroutines: 
PRT_OUTl 

Figure  24.  Subroutine  GRIDOUT  Flowchart 
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5.24  Subroutine  GT_P 
Subroutine  Call: 

GT_  P  (ii,  hrmsl,  dpi,  xktemp,  houtl,  hmax) 

Summary: 

Subroutine  GT_P  initializes  matrices  for  the  creation  of  an  internally  generated  directional 
wave  spectrum.  This  wave  spectrum  has  50  frequencies  and  36  directions. 

Input  Variables: 

dpi  Real 

ii  Integer 

hrmsl  Real 

xktemp  (points)  Real 

Output  Variables: 

hmax 
houtl 
wlen 

Local  Variables: 

hdep  Real  Breaking  Wave  Height  at  Specific  Depth 

Subroutines  Called  from  GT_P  (  ):  None. 

GT_P  (  )  Called  from  Subroutines: 

GRIDOUT 
PRT_OUTl 
PRTOUT 


Real 

Maximum  Wave  Height  / 10.0 

Real 

Significant  Wave  Height 

Real 

Wave  Length 

Offshore  Depth 

Index  where  Wave  Probabilities 
Exceed  Threshold 
Root  Mean  Square  Wave  Height 
Temporary  Variable  for  Wave  Number 
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Figure  25.  Subroutine  GT_P  Flowchart 
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5.25  Subroutine  GT_SIG_H 


Subroutine  Call: 

GT_SIG_H  (ifreq,  idirec,  esowm,  ehsig) 

Summary: 

Subroutine  GT_SIG_H  sums  the  energy  present  in  the  directional  wave  spectrum  and 
calculates  the  significant  wave  height.  The  significant  wave  height  is  defined  as: 


4  JZA/J) 


Where,  e  is  the  directional  wave  spectrum. 

Input  Variables: 

esowm  (dirNum,freqNum)  Real 
idirec  Integer 

ifreq  Integer 

Output  Variables: 

ehsig  Real 


Directional  Wave  Spectrum 

Number  of  Direction  Bins  in  Input  Spectrum 

Number  of  Frequencies  in  Input  Spectrum 


Significant  Wave  Height  from 
Directional  Spectrum 


Local  Variables: 

idir 

ifrq 

suml 

sum2 


Integer 

Integer 

Real 

Real 


Subroutines  Called  from  GT_SIG_H  (  ): 

GT_SIG_H  ( )  Called  from  Subroutines: 

CALCSURF 

READSPEC 

WAVEFTT 


Direction  Loop  Counter 
Frequency  Loop  Counter 
Summing  Variable  for  Wave  Height 
Summing  Variable  for  Wave  Height 

None. 
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Figure  26.  Subroutine  GT_SIG_H  Flowchart 


[  START  ) 


5.26  Subroutine  INITLIZE 


Subroutine  Call: 

INITLIZE  (dp,  fqd,  Cg,  xk,  c) 

Summary: 

Subroutine  INITLIZE  calculates  wave  parameters  at  the  farthest  offshore  point.  Wave  celerity 
is  calculated  from  the  dispersion  relation  given  by: 

a2-  gk  tanhffc  h) 


where,  o  is  the  angular  frequency  of  the  wave  (2jc/T),  g  is  gravity,  k  is  wave  number,  and  h  is  the  local 


Cg  =  0.5C(l+ 


2kh 
sinh  kh 


water  depth.  Wave  group  velocity  is  calculated  from  the  linear  wave  theory  relation  given  by: 
where,  C  is  the  wave  celerity. 


Input  Variables: 


dp 

fqd 


Output  Variables: 


Cg 

xk 


Real 

Real 


Real 

Real 

Real 


Offshore  Water  Depth 
Peak  Frequency 


Wave  Celerity  at  Input  Depth  &  Frequency 
Group  Velocity  at  Input  Depth  &  Frequency 
Wave  Number  at  Input  Depth  &  Frequency 
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Local  Variables: 


Real  Deep  Water  Wave  Number 

Subroutines  Called  from  INITLIZE  (  ): 

WAVENUM 

INITLIZE  (  )  Called  from  Subroutines: 

CALCS URF 

Figure  27.  Subroutine  INITLIZE  Flowchart 


[  START 


[  RETURN  ^ 


5.27  Subroutine  KLONG 


Subroutine  Call: 

KLONG  (j_ii,  xdelt_gr,  eb_last,  along,  blong,  long,  c3,  iimax,  vwind,v) 

Sununary: 

Subroutine  KLONG  calculates  longshore  current  velocity  using  an  implicit  double  sweep 
method  (Tridiagonal  Method)  based  on  the  work  of  Kraus  and  Larson  (1991).  The  central  difference 
equation  is  of  the  form: 


a,  V,-i+  bi  Vt~ a  V<+/  =  r, 


where,  V  is  the  longshore  current  velocity.  The  coefficients  a ,  b,  and  c  are  calculated  from  wave 
parameters. 

Input  Variables: 


along  (points) 

Real 

Horizontal  Mixing  Parameter 

blong  (points) 

Real 

Bottom  Friction 

c3 

Real 

Radiation  Stress  Factor  for  Longshore  Current 

clong  (points) 

Real 

Velocity 

Wind  Stress  Contribution  to 

ebjast 

Real 

Longshore  Current 

Roller  Dissipation  Term  Farthest  Offshore 

iimax 

Integer 

Number  of  Calculation  Locations 

j_ii 

Integer 

Index  where  Wave  Probabilities 

vwind 

Real 

Exceed  Threshold 

Wind  Driven  Longshore  Current  Velocity 

xdelt_gr 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

Output  Variables: 

v  (points) 

Real 

Longshore  Current  Velocity 
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Local  Variables: 


ah 

Real 

bh 

Real 

ch 

Real 

dn 

Real 

ee  (points) 

Real 

ff  (points) 

Real 

ieeff 

Integer 

ii 

Integer 

iuse 

Integer 

xde!2 

Real 

Subroutines  Called  from  KLONG  (  ): 
KLONG  ( )  Called  from  Subroutines: 


Temporary  Variable  Used  in  Longshore 

Current  Calculation 

Temporary  Variables 

Temporary  Variables 

Temporary  Variables 

Array  of  Longshore  Driving  Terms 

Array  of  Longshore  Bottom  Friction 

Array  Index 

Loop  Variable 

Array  Index  /  Loop  Variable 

Self-Adjusting  Cross-Shore  Grid  Step 

None. 


CALCSURF 


Figure  28.  Subroutine  KLONG  Flowchart 


[  START  ) 


l 


Calculate  Energy 
Nearshore 


\ 


\ 

[  RETURN  ] 
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5.28  Subroutines  LINJL 
Subroutine  Call: 

LIN_1  (ii,  dx,  dy,  x,  y) 

Summary: 

Linear  inteipolation  routine  used  to  scale  root  mean  square  wave  height  and  water  depth  to 

user-defined  grid  step  for  output  to  the  summary  text  file. 


Input  Variables: 

dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

X 

Real 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_1  ( ): 

None. 

LIN_1  (  )  Called  from  Subroutines: 


GRIDOUT 


Figure  29.  Subroutine  LIN_1  Flowchart 


(  START  ] 

I 


Linear  Interpolation  of 
Depth  or  Wave  Height 


J 

f  RETURN  ] 
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5.29  Subroutines  LIN  2 


Subroutine  Call: 


LINJ2  (ii,  dx,  dy,  x,  y) 

Summary: 

linear  interpolation  routine  used  to  scale  percent  breaking  waves  to  user-defined  grid  step  for 
output  to  the  summary  text  file. 


Input  Variables: 

dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

•  Integer 

Index  where  Wave  Probabilities 

X 

Real 

Exceed  Threshold 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_2  ( ): 

None. 

LIN_2  ( )  Called  from  Subroutines: 

■ 

GRIDOUT 


Figure  30.  Subroutine  LIN_2  Flowchart 
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530  Subroutine  LIN  3 


Subroutine  Call: 

LIN_3  (ii,  dx,  dy,  x,  y) 

Summary: 

Linear  interpolation  routine  used  to  scale  longshore  current  velocity  distribution  to  user- 
defined  grid  step  for  output  to  the  summary  text  file. 

Input  Variables: 


dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

X 

Real 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_3  ( ): 

None. 

LIN_3  (  )  Called  from  Subroutines: 


GRIDOUT 
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Figure  31.  Subroutine  LIN_3  Flowchart 
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5.31  Subroutine  LIN  4 


Subroutine  Call: 

LIN_4  (ii,  x,  y,  xi,  yi) 

Summary: 

Linear  interpolation  routine  used  to  scale  breaker  angle  to  user-defined  grid  step  for  output  to 
the  summary  text  file. 

Input  Variables: 


x  (points) 

Real 

Input  X  Value 

y  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  in  x  and  y  arrays 

xi 

Real 

Offshore  distance 

Output  Variables: 

yi 

Real 

Interpolated  Variable 

Local  Variables: 

b 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_4  ( ):  None. 


LIN_4  (  )  Called  from  Subroutines: 
GRIDOUT 
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RETURN 


532  Subroutine  LONG1 


Subroutine  Call: 

LONG1  (ii,  cl,  c2,  c3,  c4,  dp,  ebn,  hrms,  xk,  along,  blong,  clong) 
Summary: 


Subroutine  LONG1  calculates  and  outputs  longshore  current  equation  coefficients. 


Input  Variables: 

cl 

Real 

c2 

Real 

c3 

Real 

c4 

Real 

dp 

Real 

ebn 

Real 

ii 

Integer 

hrms 

Integer 

xk 

Integer 

Output  Variables: 

along  (points) 

Real 

blong  (points) 

Real 

clong  (points) 

Real 

Local  Variables: 

None. 

Eddy  Viscosity  Coefficient 

Bottom  Friction  Coefficient 

Radiation  Stress  Coefficient 

Longshore  Wind  Stress  Coefficient 

Offshore  Water  Depth 

Roller  or  Bore  Term 

Index  where  Wave  Probabilities 

Exceed  Threshold 

Root  Mean  Square  Wave  Height 

Wave  Number 


Horizontal  Mixing  Parameter 
Bottom  Friction  Parameter 
Wave  and  Wind  Parameters 


Subroutines  Called  from  LONG1  ( ):  None. 

LONG1  ( )  Called  from  Subroutines: 

MAIN_WAV 
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Figure  33.  Subroutine  LONG1  Flowchart 


f  START  ] 

1 

Calculate 

Horizontal 

Mixing 

\ 

r 

Calculate  Bottom 
Friction 

\ 

r 

Calculate  Wave 
Radiation  Stress 

1 

r 

Calculate  Wind 
Stress  Forcing 

\ 

r 

[  RETURN  ] 
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5.33  Subroutine  MAIN_WA  V 

Subroutine  Call: 

MA1N_WAV  (roller,  dxy,  xxl,  xshift,  b,  cl,  c2,  c3,  c4,  hrms,  xdelt_gr,  fqz,  nnn,  per,  xk,  fqd, 
Cg,  self_st,  dstart,  theta,  theta2,  xtemp,  xktemp,  ebjast,  htemp,  ptemp,  ebtemp,  thetatemp,  Umax] 
along,  blong,  clong,  convg,  distmax,  rk,  bl,  surf,  j_ii) 

Summary: 

Subroutine  MAIN_WAV  is  the  main  driver  for  coordinating  the  iterative  solution  method 
applied  to  solve  for  the  wave  and  current  parameters.  This  approach  is  necessary  because  several  of 
the  parameters  including  wave  height,  wave  length,  wave  celerity,  longshore  current  velocity,  and 
wave  induced  setup  are  interdependent,  as  well  as  depth  dependent. 


Input  Variables: 

b 

Real 

cl 

Real 

c2 

Real 

c3 

Real 

c4 

Real 

Cg 

Real 

dstart 

Real 

dxy  (points) 

Real 

fqd 

Real 

fqz 

Real 

hrms 

Real 

nnn 

Integer 

per 

Real 

roller 

Logical 

selfjst 

Char*l 

theta 

Real 

xdelt_gr 

Real 

xk 

Real 

xshift 

Real 

xxl  (points) 

Real 

Empirical  Factor  in  Breaking  Model  =  1.0 

Mixing/Eddy  Viscosity  Coefficient 

Bottom  Friction  Coefficient 

Factor  for  Radiation  Stress 

Friction  Coefficient  =  0.0035 

Wave  Group  Velocity 

Starting  Depth  from  Input  File 

Corresponding  Depths  with  Tide 

Peak  Frequency  at  the  Center  of  the 

Frequency  Band 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Number  of  Points  in  Input  Depth  Array 

Peak  Period  of  Directional  Wave  Spectrum 

Roller  Option  Flag  (True  or  False) 

Self  Start  Flag  (Yes  or  No) 

Wave  Angle 

Self-Adjusting  Cross-Shore  Grid  Step 
Wave  Number 

Horizontal  Cross-Shore  Location 
Adjusted  Cross-Shore  Distances  from 
Depth  Profile 
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Output  Variables: 

along  (points) 

Real 

bl  (points) 

Real 

blong  (points) 

Real 

clong  (points) 

Real 

convg 

Logical 

distmax 

Real 

eb_last 

Real 

ebtemp  (points) 

Real 

htemp  (points) 

Real 

iimax 

Integer 

j_ii 

Integer 

ptemp  (points) 

Real 

rk  (points,4) 

Real 

surf 

Logical 

theta 

Real 

theta2 

Real 

xktemp  (points) 

Real 

xtemp  (points) 

Real 

Local  Variables: 

brklO 

Logical 

cg2 

Real 

check 

Real 

conv_count 

Integer 

done 

Logical 

dp 

Real 

eb 

Real 

Horizontal  Mixing  Parameter 
Bottom  Slope 

Bottom  Friction  for  Deep  &  Shallow  Water 
Wind  Stress  Contribution  to 
Longshore  Current 
Energy  Equation  Convergence  Flag 
(True  or  False) 

Farthest  Offshore  Distance 
Roller  Dissipation  Term  at  Farthest 
Point  Offshore 

Temporary  Roller  Dissipation  Term 
Across  Transect 

Temporary  Variable  for  Significant  Wave 
Height  Values 

Number  of  Calculation  Locations 
Index  where  Wave  Probabilities 
Exceed  Threshold 
Percentage  of  Breaking  Waves  & 

Breaking  Types 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 
Flag  for  Low/No  Surf  Conditions 
(True  or  False) 

Wave  Angle 

Wave  Angle  at  Input  Starting  Depth 
Temporary  Variable  for  Wave  Number 
Temporary  Variable  for  Cross-Shore  Values 


Flag  Variable  to  Find  First  Location  Where 

10%  of  Waves  Are  Breaking  (True  or  False) 

Additional  Wave  Group  Velocity 

Difference  Between  Wave  Induced 

Setup  Calculations 

Number  of  Convergence  Iterations 

Loop  Control  Variable  for  Main  Wave 

Calculation  Loop  (True  or  False) 

Offshore  Water  Depth 
Temporary  Roller  Dissipation  Term 
Across  Transect 
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etanew  (points) 

Real 

etaold  (points) 

Real 

hrms2 

Real 

ii 

Integer 

11 

Real 

10 

Real 

P(4) 

Real 

pet (4) 

Real 

rhs 

Real 

thetaO 

Real 

xoff 

Real 

Subroutines  Called  from  MAIN_WAV  (  ): 

ABORT 

BALANCEQ 

GET_BRK 

GET_RHS 

LONG1 

PT_2 

SLF_STRT 

SETUP 

MAIN_WAV  ( )  Called  from  Subroutines: 
CALCSURF 


Wave  Induced  Setup  Estimated  at 
New  Location 

Wave  Induced  Setup  Estimated  at 
Previous  Location 

Wave  Height  for  Next  Onshore  Grid  Location 
Index  where  Wave  Probabilities 
Exceed  Threshold 

Wave  Length  at  Next  Onshore  Grid  Location 
Wave  Length  at  Grid  Cell  (1)  Offshore 
Array  for  Breaker  Percentage  Totals 
Percent  of  Different  Breaker  Types: 
pet  (1)  =  Spilling 
pet  (2)  =  Plunging 
pet  (3)  =  Surging 
pet  (4)  =  Total 

Right  Hand  Side  of  Energy  Balance  Equation 
Wave  Angle  at  Grid  (1)  Offshore 
Distance  Offshore 
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534  Subroutine  MDSRF1 
Subroutine  Call: 


MDSRF1  (alfa,  chrlie,  pet,  echo,  foxtrt,  jgamma,  ihtll,  ihtl2,  fikout) 
Summary: 


Subroutine  MDSRF1  calculates  and  prints  the  modified  surf  index  number  to  the  output  file. 


Input  Variables: 


alfa  Real 

chrlie  Real 

echo  Real 

foxtrt  Real 

ihtH  Real 

iht!2  Real 

jgamma  Integer 

Pet  (4)  Real 


file_out  Char*40 

Output  Variables:  None. 

Local  Variables: 

idir  Integer 

index  Integer 

ispd  Integer 

m  Integer 

percent  Real 

srfmod  Real 

sum  Real 

suml  Real 

sum2  Real 

temp  Real 

theta2  Real 

value  Real 


Significant  Breaker  Height 
Dominant  Breaker  Period 
Breaker  Angle 

Longshore  Current  Speed  and  Direction 
Wind  Speed 
Wind  Direction 

Temporary  Variable  set  to  Beach  Orientation 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Output  File  Name 


Index  for  Surf  Index  Wind  Direction 
Breaker  Type  Indicator  for  Surf  Index 
Index  for  Surf  Index  Wind  Speed  Lookup  in 
Modification  Table 

Temporary  Variable  to  Rotate  Direction 
Breaker  Type  Percentage 
Modified  Surf  Index  from  Sum  of  Values 
Resulting  from  Navy  Modification  Tables  in 
MDSRF2  ( ) 

Running  Total  of  Surf  Index 
Modified  Surf  Index  Value  for  Wave  Angle 
Value  for  Longshore  Current 
Temporary  Wave  Angle  Variable 
Rotated  Wind  Direction 
Modification  Number 
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wind  (3,3,8)  Real  Surf  Index  Wind  Modification  Table 

Subroutines  Called  from  MDSRF1  (  ): 

MDSRF2 

MDSRF1  ( )  Called  from  Subroutines: 

SURF 
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Figure  35.  Subroutine  MDSRF1  Flowchart 
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5.35  Subroutine  MDSRF2 


Subroutine  Call: 

MDSRF2  (index,  xin,  yin,  value) 

Summary: 

Subroutine  MDSRF2  contains  the  modified  surf  index  (MSI)  tables.  The  MSI  number 
calculated  using  a  two  dimensional  linear  interpolation  by  areas. 


Input  Variables: 

index 

Integer 

xin 

Real 

yin 

Real 

Output  Variables: 

value 

Real 

Local  Variables: 

i 

Integer 

il 

Integer 

i2 

Integer 

ii 

Integer 

ix(4) 

Real 

jy(4) 

Real 

j 

Integer 

jl 

Integer 

J2 

Integer 

jj 

Integer 

tempi 

Real 

x(ll) 

Real 

xO  (4,1 1) 

Real 

xdata 

Real 

y(H) 

Real 

y0(4,ii) 

Real 

ydata 

Real 

z  (11,11) 

Real 

zO  (4,11,11) 

Real 

zl  (40)-zll(40) 

Real 

Indicator  of  Breaker  Type 
X-Coordination  for  Surf  Index 
Modification  Matrix 
Y-Coordination  for  Surf  Index 
Modification  Matrix 


Returns  Modified  Surf  Index  Number 


Loop  Counter  or  Array  Index 

Loop  Counter  or  Array  Index 

Loop  Counter  or  Array  Index 

Loop  Counter  or  Array  Index 

All  Values  Set  to  1 1.00 

Values  Set  to  10.0, 11.0, 11.0, 9.0 

Loop  Counter  or  Array  Index 

Loop  Counter  or  Array  Indent 

Loop  Counter  or  Array  Index 

Loop  Counter  or  Array  Index 

Temporary  Variable  Used  for  Interpolation 

MSI  Indices 

Breaker  Period  Modification  table 
Temporary  Index 
MSI  Indices 

Wave  Angle  Modification  table 
Temporary  Index 
Breaker  Modification  Matrix 
Whole  Breaker  Modification  Matrix 
Partial  Breaker  Modification  Arrays 


zl2  (44) 
zzO  (484) 


Real 

Real 


Partial  Breaker  Modification  Array 
Equivalent  to  zO  (1,1,1)  zzO  (1) 


Subroutines  Called  from  MDSRF2  ( ):  None. 

MDSRF2  ( )  Called  from  Subroutines: 

MDSRF1 

Figure  36.  Subroutine  MDSRF2  Flowchart 
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5.36  Subroutine  NEW_BRK 
Subroutine  Call: 

NEW_BRK  (Umax,  bl,  rk,  htemp,  wid_ii,  p2) 

Summary: 

Subroutine  NEW_BRK  calculates  a  new  percentage  of  breaker  types  from  the  highest  10%  of 
the  wave  heights  (hrms)  when  the  bottom  slope  is  positive. 


Input  Variables: 


bl  (points) 

Real 

Bottom  Slope 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

Umax 

Integer 

Number  of  Calculation  Locations 

rk  (points,4) 

Real 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 

wid_U 

Integer 

Offshore  Location  for  Surf  Zone  Width 

Output  Variables: 

p2  (4) 

Real 

Percent  of  Different  Breaker  Types  - 

Equivalent  to  pet  (4) 

p2(l)  =  Spilling 

p2  (2)  =  Plunging 

p2  (3)  =  Surging 

p2  (4)  =  Total 

Local  Variables: 

akl  (points) 

Real 

Temporary  Array  for  Wave  Height 

bkl  (points) 

Real 

Temporary  Array  Breaker  Type  =  1  SpilUng 

bk2  (points) 

Real 

Temporary  Array  Breaker  Type  =  2  Plunging 

bk3  (points) 

Real 

Temporary  Array  Breaker  Type  =  3  Surging 

bk4  (points) 

Real 

Temporary  Array  Breaker  Type  =  4  Total 
Total  Percentage  of  Breakers 

i 

Integer 

Loop  Counter 

ii 

Integer 

Loop  Counter 

nval 

Integer 

Number  of  Positive  Slope  Occurrences 

xl 

Real 

0. 1  %  of  Highest  Breakers  to  Examine  for  Type 
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x2 


Integer 


Loop  Limit  -  Set  to  Top  Percentage  of 
Significant  Wave  Height  Values 


Subroutines  Called  from  NEWJBRK  ( ): 

ABORT 

RN2 

NEW_BRK  (  )  Called  from  Subroutines: 
SHORTOUT 
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Figure  37.  Subroutine  NEW_BRK  Flowchart 
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537  Subroutine  PERCENT 


Subroutine  Call: 

PERCENT  (hrms,  period,  dp,  slope,  p) 

Summary: 

Subroutine  PERCENT  calculates  the  percentage  of  each  type  of  breaking  wave  in  the  surf 


zone. 


Input  Variables: 


dp 

Real 

hrms 

Real 

period 

Real 

slope 

Real 

Output  Variables: 

P(4) 

Real 

Local  Variables: 

firac  (3) 

Real 

gtemp 

Real 

hhigh 

Real 

hlow 

Real 

integrat 

Real 

P_flag 

Logical 

param 

Real 

Subroutines  Called 

from  PERCENT  ( ): 

GET  P 

Offshore  Water  Depth 
Root  Mean  Square  Wave  Height 
Peak  Period 
Bottom  Slope 


Array  of  Percentage  of  Breaker  Types 

pet  (1)  -  Spilling 

pet  (2)  -  Plunging 

pet  (3)  -  Surging 

pet  (4)  -  Total  Percentage 


Array  for  Percentage  Breaker  Totals 
Gravity 

Upper  Bound  of  Integration 
Lower  Bound  of  Integration 
Wave  Height  Distribution  Calculated  at 
Specific  Location 

Weighting  Factor  Flag  (True  or  False) 
Integral  Multiplier 
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Functions  Called  from  PERCENT  (  ): 
INTEGRAT 


PERCENT  (  )  Called  from  Subroutines: 

GETJBRK 

SLF_STRT 

Figure  38.  Subroutine  PERCENT  Flowchart 


f  START  ) 

M 

<  Setting 
Integral 
Multiplier 


I 


[  RETURN  ] 
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5.38  Subroutine  PRT_OUTl 
Subroutine  Call: 


PRT.OUT1  (j_ii,  xdelt,  iimax,  dxy,  xtemp,  xktemp,  thetatemp,  htemp,  ptemp,  v) 
Summary: 

Subroutine  PRT_OUTl  prints  columnar  data,  cross-shore  distributions  of  wave  and  surf 
parameters,  to  the  detailed  SURF  output  file  when  requested  by  the  user.  This  data  is  interpolated  to 


the  user  defined  grid  step,  if  possible. 

Input  Variables: 


dxy  (points)  Real 

J-1*  Integer 

iimax  Integer 

htemp  (points)  Real 

ptemp  (points)  Real 

v  (points)  Real 

xdelt  Real 

xktemp  (points)  Real 

xtemp  (points)  Real 

Output  Variables:  None. 

Local  Variables: 

dp!  Real 

hmax  Real 

houtl  Real 

hrmsl  Real 

“  Integer 

II  Integer 

pbreak  Real 

vlngl  Real 

wlen  Real 

xoffl  Real 


Corresponding  Depths  with  Tide 
Index  where  Wave  Probabilities 
Exceed  Threshold 
Number  of  Calculation  Locations 
Temporary  Variable  for  Significant  Wave 
Height  Values 

Percentage  of  Breaking  Waves  & 

Breaker  Types 

Longshore  Current  Velocity 

Surf  Zone  Output  Interval 

Temporary  Wave  Number  Array 

Temporary  Variable  for  Cross-Shore  Values 


Offshore  Depth 
Maximum  Wave  Height 
Significant  Wave  Height 
Root  Mean  Square  Wave  Height 
Array  Index  Number 
Iteration  Count 
Percentage  Breaking  Waves 
Longshore  Current  Velocity 
Wave  Length 
Distance  Offshore 
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Subroutines  Called  from  PRT_OUTl  (  ): 

ABORT 

GT_P 

GRIDOUT 

PRT_OUTl  ( )  Called  from  Subroutines: 
CALCSURF 

Figure  39.  Subroutine  PRT_OUTl  Flowchart 


i 
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539  Subroutine  PRT_OUT2 
Subroutine  Call: 


1 0UT2  (j_ii,  xdelt,  iimax,  dxy,  xtemp,  xktemp,  thetatemp,  htemp,  ptemp,  v) 

Summary: 

Subroutine  PRT_OUT2  writes  the  detailed  surf  output. 


Input  Variables: 


dxy  (points)  Real 

j_ii  Integer 

iimax  Integer 

htemp  (points)  Real 

ptemp  (points)  Real 

v  (points)  Real 

xdelt  Real 

xktemp  (points)  Real 

xtemp  (points)  Real 

Output  Variables:  None. 

Local  Variables: 

dpi  Real 

hmax  Real 

houtl  Real 

hnnsl  Real 

ii  Integer 

jj  Integer 

pbreak  Real 

vlng  Real 

wlen  Real 

xoffl  Real 


Subroutines  Called  from  PRT_OUT2  ( ): 

GT_P 

ABORT 


Corresponding  Depths  with  Tide 
Index  where  Wave  Probabilities 
Exceed  Threshold 
Number  of  Calculation  Locations 
Temporary  Variable  for  Significant  Wave 
Height  Values 

Percentage  of  Breaking  Waves  and 

Breaker  Types 

Longshore  Current  Velocity 

Surf  Zone  Output  Interval 

Temporary  Wave  Number  Array 

Temporary  Variable  for  Cross-Shore  Values 


Offshore  Depth 
Maximum  Wave  Height 
Significant  Wave  Height 
Root  Mean  Square  Wave  Height 
Array  Index  Number 
Iteration  Counter 
Percentage  Breaking  Waves 
Longshore  Current  Velocity 
Wave  Length 
Distance  Offshore 
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PRT_OUT2  (  )  Called  from  Subroutines: 
CALCSURF 

Figure  40.  Subroutine  PRT_OUT2  Flowchart 


Do  While  Offshore  Depth  is  > 
1.0  &  While  Loop  Count  is  < 
Number  of  Grid  Points 


\ 


Convert  Units  of  Distance 
Offshore,  Wave  Height  & 

I  Longshore  Current  Velocity  I 


Set  Percentage 
Breaking  Wave 


I 


Write  Detailed  j 
SURF  Output 


\ 


Calculate  Next  Grid  Point 
Distance  Offshore 


(  RETURN  ] 


101 


5.40  Subroutine  PRT_OUT3 
Subroutine  Call: 


PRT_OUT3  (file_dat) 

Summary: 

Subroutine  PRT_OUT3  writes  out  the  detailed  output  from  the  model. 

Input  Variables: 

flle-dat  Char*40  Output  File  name  *.dat 

Output  Variables:  None. 

Local  Variables: 

^ne  Char*80  Temporary  String 

Subroutines  Called  from  PRT_OUT3  ( ):  None. 

PRT_OUT3  ( )  Called  from  Subroutines: 

SURF 

Figure  41.  Subroutine  PRT_OUT3  Flowchart 


C 


START 


I 


J 


Writing 

Detailed 

Output 


5.41  Subroutine  PT2 


Subroutine  Call: 

PT2  (10,  thetaO,  fqd,  dp,  theta,  xk,  1,  Cg) 

Summary: 

Subroutine  PT2  calculates  wave  parameters  from  linear  theory  relations. 


Cg  =  nC 

Group  Velocity 

n  =  -  [1  + 

2  kh 

2 

sinh  2  kh 

sin  0 

sin  0o 

C 

Co 

Wave  angle  from  Snell’s  law 

Input  Variables: 

dp 

Real 

Offshore  Water  Depth 

fqd 

Real 

Peak  Frequency 

10 

Real 

Wave  Length  at  Offshore  Point 

thetaO 

Real 

Wave  Angle  at  Offshore  Point 

xk 

Real 

Wave  Number 

Output  Variables: 

Cg 

Real 

Group  Velocity 

1 

Real 

Wave  Length 

theta 

Real 

Wave  Angle 

xk 

Real 

Wave  Number 

Local  Variables: 

c 

Real 

Temporary  Variable 

Subroutines  Called  from  PT2  ( ): 


WAVENUM 

PT2  ( )  Called  from  Subroutines: 
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MAIN_WAV 

SLFJSTRT 

Figure  42.  Subroutine  PT2  Flowchart 


£  START  ] 

I 


Call 

WAVENUM 


I 


Calculate  Wave 
Parameters 


I 

[  RETURN  ] 
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5.42  Subroutine  RADJST1 


Subroutine  Call: 

RAD_ST1  (ifreq,  freq,  idirec,  xfrom,  esowm,  freql,  freq2,  dstart,  igamma,  theta,  hrms,  surf, 
fqd,  per,  fqz) 

Summary: 

Subroutine  RADJST1  searches  the  directional  wave  spectrum  to  identify  the  dominant  wave 
frequency  and  sums  the  wave  energy  directed  toward  shore.  The  flux  of  momentum  or  Radiation 
Stress,  which  contributes  to  driving  the  longshore  current,  is  calculated  following  Thornton  and  Guza 

Sxy(O.'f)  =  E ( 0,  f) n (f)  sin  a(f)  cos  a(f) 


(1986). 

In  the  above  equation  Sxy  is  the  Radiation  Stress,  E  is  the  total  energy  in  the  directional  wave 
spectrum,  n  is  the  ratio  of  wave  group  velocity  to  wave  velocity,  and  a  is  the  wave  angle.  The  ratio  n 
from  linear  wave  theory  is  given  by: 


n  = 


Q 

C 


=  0.5(1+ 


2kh 
sinh  kh 


where,  C%  is  the  group  velocity,  C  is  the  wave  velocity  or  celerity,  k  is  the  wave  number  and  h  is  the 
local  water  depth. 

Input  Variables: 


dstart 

Real 

Input  Starting  Depth 

esowm  (dirNum,freqNum) 

Real 

Directional  Wave  Spectrum 

freq  (fteqNum) 

Real 

Input  Wave  Spectrum  Center  Frequencies 

freql(freqNum) 

Real 

Beginning  Frequency  Bin  Values 

freq2  (freqNum) 

Real 

Ending  Frequency  Bin  Values 

idirec 

Integer 

Number  of  Directions  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequencies  in  Input  Spectrum 

igamma 

Integer 

Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 

xfrom  (freqNum) 

Real 

Direction  Array,  Direction  Wave  Energy 
Comes  From 
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Output  Variables: 

fqd 

Real 

fqz 

Real 

hrms 

Real 

per 

Real 

surf 

Logical 

theta 

Real 

Local  Variables: 

direc 

Real 

ees 

Real 

esum 

Real 

esumm 

Real 

ffd 

Real 

idir 

Integer 

ifrq 

Integer 

m 

Integer 

maxfrq 

Integer 

summax 

Real 

sumzero 

Real 

sxy 

Real 

sxysum 

F 

temp 

Real 

temp2 

Real 

theta2 

Real 

xk 

Real 

xkd 

Real 

Subroutines  Called  from  RADJST1  ( ): 

RADJST2 

WAVENUM 

RAD_ST1  ( )  Called  from  Subroutines: 
CALCSURF 


Peak  Frequency  at  the  Center  of  the 

Frequency  Band 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 

Flag  for  Low  or  No  Surf  Conditions 

(True  or  False) 

Wave  Angle 


Wave  Direction 

Spectral  Density  at  a  Particular 

Frequency  and  Direction 

Sum  of  Energy  in  One  Frequency  Band 

Over  all  Directions 

Sum  of  All  Energy  in  Directional  Spectrum 
Wave  Frequency 
Loop  Counter 
Loop  Counter 

Temporary  Variable  for  Rotating  Wave  Angle 
Frequency  at  Maximum  Spectral  Density 
Frequency  Band  with  Maximum  Energy 
Summation  of  Zero-Crossing  Frequency 
Energy 

Radiation  Stress 

Sum  of  Radiation  Stress  Energy 
Temporary  Variable  in  Radiation 
Stress  Calculation 

Temporary  Variable  for  Frequency  Band  with 
Maximum  Energy 

Angle  Between  Wave  Ray  and  Beach 
Perpendicular  Projection 
Wave  Number 

Wave  Number  Multiplied  by  the 
Local  Water  Depth 
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5.43  Subroutine  RADST2 
Subroutine  Call: 


RAD_ST2  (freq,  sxysum,  sumzero,  esumm,  maxfrq,  dstart,  theta,  hrms,  surf,  fqd,  per,  fqz) 

Summary: 

Subroutine  RAD_ST2  calculates  several  parameters  based  on  the  total  energy  in  the 
directional  wave  spectrum.  A  check  is  performed  to  confirm  that  wave  energy  is  directed  onshore 
before  writing  summary  information  to  the  output  file. 

Input  Variables: 


dstart 

esumm 

freq  (freqNum) 

maxfrq 

sumzero 

sxysum 

Output  Variables: 


Real 

Real 

Real 

Integer 

Real 

Real 


Input  Starting  Depth 

Sum  of  All  Energy  in  Directional  Spectrum 
Input  Wave  Spectrum  Center  Frequencies 
Frequency  at  Maximum  Spectral  Density 
Summation  of  Zero-Crossing 
Frequency  Energy 

Sum  of  Radiation  Stress  energy 


fqd 

fqz 

hrms 

per 

surf 

theta 

Local  Variables: 


Real 

Real 

Real 

Real 

Logical 

Real 


Peak  Frequency 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 

Logical  Flag  for  Low/No  Surf  Conditions 

(True  or  False) 

Wave  Angle 


hs 

sxy2 

temp 

theta3 

xk 

xkd 


Real 

Real 

Real 

Real 

Real 

Real 


Significant  Wave  Height 

Temporary  Wave  Energy 

Temporary  Variable  for  Energy 

Wave  Angle  in  Degrees 

Wave  Number  Calculated  at  Peak  Frequency 

and  Input  Starting  Depth 

Wave  Number  *  Water  Depth 


Subroutines  Called  from  RAD_ST2  (  ): 

S_NOSURF 

WAVENUM 

RADJST2  (  )  Called  from  Subroutines: 
RAD.ST1 

Figure  44.  Subroutine  RAD_ST2  Flowchart 
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5.44  Subroutine  READRFRC 
Subroutine  Call: 


READRFRC  (fracname,  ifreq,  ffeq,  idirec,  xfrom,  xcoeff,  xtheta) 

Summary: 

Subroutine  READRFRC  reads  refraction  information  from  a  formatted  input  file.  The 
matrices  contained  in  these  files  are  used  to  shoal  and  refract  a  directional  wave  spectrum  from  an 
offshore  point  to  a  location  where  depth  information  is  available.  The  number  of  frequency  bins  must 
notexceed  50  and  the  number  of  direction  bins  must  not  exceed  180.  The  directional  coverage  of  the 
refraction  and  shoaling  coefficients  must  range  from  0  to  360  degrees.  Partial  coverage  over  a 
fraction  of  the  compass  (e.g.  180  degree  sector)  will  introduce  errors. 


Input  Variables: 

fracname 

Char*40 

Output  Variables: 

idwsdirec 

Integer 

idwsfreq 

Integer 

sdir  (dirNum) 

Real 

sfreq  (freqNum) 

Real 

xcoeff  (dirNum,freqNum) 

Real 

xtheta  (dirNum,freqNum) 

Real 

Local  Variables: 

cfmatch 

Logical 

cfreq  (freqNum) 

Real 

col 

Real 

dangle 

Real 

dir 

Real 

dirin 

Integer 

dirord 

Integer 

dirouts  (dirNum) 

Real 

dirs  (dirNum) 

Real 

dmatch 

Logical 

Wave  Refraction  File 


Number  of  rows  (Directions)  in  the 
Directional  Wave  Spectrum  Matrix 
Number  of  columns  (Frequencies)  in  the 
Directional  Wave  Spectrum  Matrix 
Direction  Array  for  each  bin  in  the 
Directional  Wave  Spectrum 
Center  Frequency  of  each  Directional 
Wave  Spectrum 

Wave  Height  Refraction  Coefficients 
Angle  Refraction  Coefficients 


Flag  for  Center  Frequency  Match 
Center  Frequency  of  each  Bin 
Number  of  Columns 
Angle  Between  Directional  Bins 
Number  of  Angles 
X-Coordinates  of  known  values 
Direction  of  Waves 

1  -  Direction  Waves  are  coming  from 

2  -  Direction  Waves  are  going  to 
Interpolated  X-Coordinates 

Temporary  Direction  Wave  Energy  Comes  Rom 
Rag  for  Directional  Match 


ItO 


dots 

Integer 

drl 

Real 

dth 

Real 

dum 

Real 

dum2 

Real 

dumstr 

Char*80 

fmatch 

Logical 

fnum 

Integer 

found 

Integer 

frchk 

Integer 

frq 

Real 

I 

Integer 

ii 

Integer 

icol 

Integer 

idir 

Integer 

idirec 

Integer 

ifreq 

Integer 

ifrq 

Integer 

instat 

Integer 

irow 

Integer 

j 

Integer 

jj 

Integer 

k 

Integer 

kk 

Integer 

Ifreq 

Real 

lowcut 

Integer 

mpnt 

Integer 

refs  (dirNum) 

Real 

rfrtmp  (dirNum,freqNum) 

Real 

row 

Real 

rtmpout  (dirNum) 

Real 

sfreqin  (dirNum) 

Real 

shltmp  (dirNum, freqNum) 

Real 

splout  (dirNum) 

Real 

stmpout  (dirNum) 

Real 

temp  (dirNum,freqNum) 

Real 

temp2  (dirNum,ffeqNum) 

Real 

tmpinr  (dirNum) 

Real 

tmpins  (dirNum) 

Real 

ufreq 

Real 

upcut 

Integer 

xfrom  (dirNum) 

Real 

Y-Coordinates  of  known  values 
Initial  Direction  Bin 

Temporary  Angle  Between  Directional  Bins 

Temporary  Variable 

Temporary  Variable 

Temporary  Variable 

Flag  for  Frequency  Match 

Bin  Number 

Flag  Indicator 

Total  Number  of  Frequencies 
Number  of  Frequencies 
Loop  Counter 
Counter 

Number  of  Columns 
Loop  Counter 

Number  of  Rows  (Directions)  in  the 

Refraction/Shoaling  Matrix 

Number  of  Columns  (Frequencies)  in  the 

Refraction/Shoaling  Matrix 

Loop  Counter 

Error  Status 

Number  of  Rows 

Loop  Counter 

Counter 

Counter 

Counter 

Lower  Frequency  Bin  Limit 
Lower  Cut  Off  Limit 
Number  of  Rows  divided  by  2 
Temporary  Array 

Temporary  Matrix  for  Reversing  Wave 

Direction 

Number  of  Rows 

Interpolated  Coordinates 

Temporary  Frequency  Array 

Temporary  Matrix  for  Reversing  Wave 

Direction 

interpolated  Y-Coordinates 

Interpolated  Coordinates 

Temporary  Variable 

Temporary  Variable 

Temporary  Variable 

Temporary  Variable 

Upper  Frequency  Bin  Limit 

Upper  Cut  Off  Limit 

Direction  Wave  Energy  Comes  From 


lit 


Subroutines  Called  from  READRFRC  (  ): 

ABORT 

GENRLSPL 

READRFRC  ( )  Called  from  Subroutines: 
SURF 
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Figure  45.  Subroutine  READRFRC  Flowchart 


5.45  Subroutine  READSPEC 


Subroutine  Call: 


READSPEC  (ifreq,  idirec,  Cfreq,  Lfreq,  Ufreq,  xfrom,  sowm,  period,  ehsig,  dangle,  spefile) 

Summary: 


Subroutine  READSPEC  opens  and  reads  a  directional  wave  spectrum  file,  which  must 
conform  to  a  specific  format,  but  the  number  of  frequencies  and  directions  can  vary.  The  maximum 
number  of  directions  is  180  and  the  maximum  number  of  frequencies  is  50.  The  directions  should  be 
evenly  spaced,  and  the  frequency  bins  can  be  fixed  or  variable  width  with  units  of  energy  density 
(mA2/(Hz*radians)).  This  energy  density  matrix  is  initialized,  filled,  and  converted  to  units  of  feet 


squared  inside  this  subroutine.  In  addition,  the  direction  of  wave  energy  can  be  the  direction  FROM 
which  waves  are  coming  or  TO  which  waves  are  going  as  denoted  in  the  tenth  header  line  by  a  1  or  2 
respectively.  The  directional  wave  spectrum  must  be  defined  from  0  to  360  degrees.  Use  of  partial 
directional  sectors  (e.g.  0  to  180  degrees)  will  cause  errors. 


Input  Variables: 

None. 

Output  Variables: 

Cfreq  (fireqNum) 

Real 

dangle 

Real 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

idirec 

Integer 

ifreq 

Integer 

Lfreq  (freqNum) 

Real 

period  (freqNum) 

Real 

spefile 

Char*40 

Ufreq  (freqNum) 

Real 

xfrom  (dirNum) 

Real 

Center  Frequency  Bin  Limit 

Angle  Between  Directional  Bins 

Significant  Wave  Height  from 

Directional  Spectrum 

Directional  Wave  Spectrum 

Number  of  Direction  Bins  in  Input  Spectrum 

Number  of  Frequency  Bins  in  Input  Spectrum 

Lower  Frequency  Bin  Limit 

Period  Array  (1 /Frequency) 

Wave  Spectrum  File  Name 
Upper  Frequency  Bin  Limit 
Direction  Array,  Direction  Wave  Energy 
Comes  From 


114 


Local  Variables: 


col 

Real 

df 

Real 

dir 

Real 

dirord 

Integer 

dth 

Real 

dum 

Char*l 

drl 

Real 

fnum 

Integer 

frq 

Real 

ftsq2msg 

Real 

I 

Integer 

icol 

Integer 

idir 

Integer 

ifrq 

Integer 

instat 

Integer 

irow 

Integer 

j 

Integer 

mpnt 

Integer 

mult 

Real 

row 

Real 

temp  (dirNum,dirNum) 

Real 

Number  of  Columns 

Difference  between  Upper  &  Lower  Bins 
Number  of  Angles 
Direction  of  Waves 

1  -  Direction  Waves  are  coming  from 

2  -  Direction  Waves  are  going  to 
Width  of  Direction  Bin 
Temporary  Variable 

Initial  Direction  Bin 
Bin  Number 
Number  of  Frequencies 
Conversion  Factor 
Loop  Counter 
Number  of  Columns 
Direction  Loop  Counter 
Loop  Counter 
Error  Status 
Number  of  Rows 
Loop  Counter 

Number  of  Rows  divided  by  2 
Temporary  Calculation  Variable 
Number  of  Rows 
Temporary  Array 


Subroutines  Called  from  READSPEC  ( ): 

ABORT 

GT_SIG_H 

READSPEC  (  )  Called  from  Subroutines: 
SURF 
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Figure  46.  Subroutine  READSPEC  Flowchart 


116 


5.46  Subroutine  REFRAC 


Subroutine  Call: 

REFRAC  (idirec,  ifreq,  xffom,  xtheta,  xcoeff,  esowm,  ehsig) 

Summary: 

For  each  frequency  and  direction  bin  in  the  input  directional  wave  spectrum,  the  shallow  water 
direction  band  for  each  deep  water  direction  band  is  found.  Wave  energy  from  each  deep  water  band 
is  multiplied  by  the  combined  reffaction/shoaling  coefficient  and  moved  into  the  proper  shallow  water 
band  to  provide  a  shallow  water  directional  spectrum. 


Input  Variables: 


dangle 

Real 

Angle  Between  Directional  Bins 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Real 

Number  of  Frequencies  in  Input  Spectrum 

xcoeff  (dirNum,ffeqNum) 

Real 

Wave  Height  Refraction  Coefficients 

xtheta  (dirNum,ffeqNum) 

Real 

Angle  Refraction  Coefficients 

Output  Variables: 

ehsig 

Real 

Significant  Wave  Height  from 

Directional  Spectrum 

esowm  (dirNum,freqNum) 

Real 

Directional  Wave  Spectrum 

Local  Variables: 

esite  (dirNum,freqNum) 

Real 

Directional  Spectrum  in  Shallow  Water 

idir 

Integer 

Direction  Loop  Counter 

ifrq 

Integer 

Frequency  Loop  Counter 

itemp 

Integer 

Temporary  Wave  Angle  Variable 

itheta  (dirNum,ffeqNum) 

Integer 

Shoreward  Energy  Spectrum 

jdir 

Integer 

Loop  Variable 

mtemp 

Integer 

Temporary  Wave  Angle  Variable 

sum 

Real 

Temporary  Wave  Energy  Summation  Variable 

sum2 

Real 

Temporary  Wave  Energy  Summation  Variable 

ytheta 

Real 

Temporary  Wave  Angle  Variable 

Subroutines  Called  from  REFRAC  ( ):  None. 
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REFRAC  (  )  Called  from  Subroutines: 
SURF 

Figure  47.  Subroutine  REFRAC  Flowchart 
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Direction  Index  for  All 
Deep  Water  Components 


H  End  of  Frequency  Loop 
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I 


Number  of 

_ |Direction  Bins  in| 

Input  Spectrum 
Loop 


Number  of 
jFrequency  Bands 
in  Input 
Spectrum  Loop 


Set  Sum,  Initialize  to  a 
Nonzero  Value 


i 


Number  of 
-^Direction  Bins  in 
Inner  Loop 


j 


Sum  Energy  in 
Each  Direction 
Bin 


\ 


End  of  Loop 


I 


Number  of 
Frequency  Bands 
Loop 


\ 


Number  of 
mmi  Direction  Bins 
Loop 


Calculate  Frequency 
Spectrum 


[ 

[  RETURN  ] 
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5.47  Subroutine  RN2 


Subroutine  Call: 

RN2  (n,  x,  yl,  y2,  y3,  y4) 

Summary: 

Subroutine  RN2  calculates  percentages  of  each  type  of  breaker  in  the  surf  zone. 


Input  Variables: 


n 

Integer 

Number  of  Waves  Considered  Breaking  on  a 
Positive  Bottom  Slope 

x  (points) 

Real 

Temporary  Significant  Wave  Height  Array 

yl  (points) 

Real 

Spilling  Breaker  Type 

y2  (points) 

Real 

Plunging  Breaker  Type 

y3  (points) 

Real 

Surging  Breaker  Type 

y4  (points) 

Real 

Total  Number  of  Breakers 

Output  Variables: 

yl  (points) 

Real 

Spilling  Array  Breaker  Type 

y2  (points) 

Real 

Plunging  Array  Breaker  Type 

y3  (points) 

Real 

Surging  Array  Breaker  Type 

y4  (points) 

Real 

Total  Array  Breaker  Type 

Local  Variables: 

hold 

Real 

Temporary  Variable  Used  for  Repositioning 

i 

Integer 

Loop  Counter 

j 

Integer 

Loop  Counter 

js 

Integer 

Loop  Starting  Index 

m 

Integer 

Number  of  Waves  Considered  Breaking  on  a 
Positive  Slope 

Subroutines  Called  from  RN2  (  ): 

None. 

RN2  ( )  Called  from  Subroutines: 
NEWJBRK 
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Figure  48.  Subroutine  RN2  Flowchart 


[  START  ] 


I 


5.48  Subroutine  S  COEFF 


Subroutine  Call: 

S_COEFF  (dp,  fqd,  hrms,  theta,  c,  xk,  wdir,  igamma,  wdspd,  cl,c2,  c3,  c4,  cf,  vwind) 

Summary: 

Subroutine  S_COEFF  calculates  several  parameters  in  the  longshore  current  equation 
including  the  Radiation  Stress,  the  bottom  stress,  and  the  wind  stress.  A  check  is  performed  to  assure 
that  wave  induced  motion  is  not  dominated  by  wind  effects  and  a  warning  message  is  written  to  the 
output  file  if  this  condition  is  violated.  An  assumption  is  made  that  if  the  wave  induced  orbital 
velocity  is  greater  than  the  wind-forced  component  of  the  longshore  current,  the  local  conditions  are 
wave  dominated. 


Input  Variables: 


c 

Real 

Wave  Celerity  at  Input  Starting  Depth 

dp 

Real 

Water  Depth  Offshore 

fqd 

Real 

Peak  Frequency  from  Directional  Spectrum 

hrms 

Real 

Root  Mean  Square  Wave  Height 

igamma 

Integer 

Beach  Orientation,  Compass  Heading 
Directly  Toward  Beach 

theta 

Real 

Wave  Angle 

wdir 

Real 

Input  Wind  Direction  Compass  Heading 

wdspd 

Real 

Input  Wind  Speed 

xk 

Real 

Wave  Length  at  Input  Starting  Depth 

Output  Variables: 

cl 

Real 

Mixing/Eddy  Viscosity  Coefficient 

c2 

Real 

Bottom  Friction  Coefficient 

c3 

Real 

Factor  for  Radiation  Stress 

c4 

Real 

Friction  Coefficient 

vwind 

Real 

Wind  Driven  Longshore  Current  Velocity 
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Local  Variables: 


c4tmp 

cd 

cf 

dwind 

m 


theta4 

uoib 

xn 


Real 

Real 

Real 

Real 

Integer 


Real 

Real 

Real 


Subroutines  Called  from  S_COEFF  (  ): 
S_COEFF  (  )  Called  from  Subroutines: 
CALCSURF 


Temporary  Variable  Used  in  Wind  Velocity 

Vector  Calculation 

Coefficient  of  Drag  Used  in  Wind 

Velocity  Calculation 

Coefficient  of  Friction  for  the  Bottom 

Sign  of  Wind  Vector  (Positive  or  Negative) 

Temporary  Variable  Used  in  Rotating 

Wind  Angle 

Rotated  Wind  Direction 

Wave  Particle  Orbital  Velocity 

Eddy  Viscosity  Mixing  Coefficient 

None. 
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Figure  49.  Subroutine  S_COEFF  Flowchart 


[  START  ] 
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5.49  Subroutine  S.NOSURF 
Subroutine  Call: 

S_NOSURF  (  hsig,  surf  ) 

Summary: 

Subroutine  S_NOSURF  is  called  to  determine  if  local  conditions  are  significant  enough  to 
proceed  with  surf  zone  calculations.  The  minimum  condition  for  continuation  is  that  the  significant 
wave  height  calculated  from  the  directional  wave  spectrum  must  be  greater  than  0.15  m. 

Significant  Wave  Height 

Flag  to  Indicate  Low  or  No  Surf  Conditions 
(True  or  False) 

Subroutines  Called  from  S_NOSURF  ( ):  None. 

S_NOSURF  ( )  Called  from  Subroutines: 

CALCSURF 

RADJST2 


Input  Variables: 

hsig  Real 

Output  Variables: 

surf  Logical 

Local  Variables:  None. 
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Figure  50.  Subroutine  S_NOSURF  Flowchart 


Set  Surf  Flag 
to  True 


Set  Surf  Flag 
to  False 


5.50  Subroutine  S_TIDE 


Subroutine  Call: 

S_TIDE  (tide,  ydepth,  nnn,  dxyl,  xxl,  dxy,  xshift) 

Summary: 

Subroutine  S_TIDE  adds  the  tidal  elevation  to  each  cross-shore  point  in  the  input  depth 

profile. 


Input  Variables: 


dxyl  (points)  Real 

nnn  Integer 

tide  Real 

xxl  (points)  Real 

ydepth  Char*l 

Output  Variables: 

dxy  (points)  Real 

xshift  Real 

Local  Variables: 

ddiff  Real 

n  Integer 

nn  Integer 

111111  Integer 

xdiff  Real 

ztide  Real 


Subroutines  Called  from  S_TIDE  (  ): 
S_TIDE  (  )  Called  from  Subroutines: 
CALCS URF 


Corresponding  Depths  without  Tide 
Number  of  Points  in  Input  Depth  Array 
Tide  Level 

Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

Usage  of  Input  Depth  (Yes/No) 


Adjusted  Depths  with  Tide 
Offshore  Distance 


Change  in  Water  Depth 
Loop  Counter 
Loop  Counter 
Loop  Counter 

Change  in  Cross-Shore  Location 
Tide  Level 

None. 
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Figure  51.  Subroutine  S_TIDE  Flowchart 
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5.51  Subroutine  SEAFIT 


Subroutine  Call: 

SEAFIT  (hsig,  per,  dir,  ifreq,  idirec,  freql,  freq2,  xfrom,  esowm) 

Summary: 

Subroutine  SEAFIT  calculates  a  directional  wave  spectrum  from  an  input  wave  height  and 
wave  period  using  a  Pierson-Moskowitz  spectrum  representation  and  a  cosine  to  the  fourth  directional 
spreading  function.  The  modified  Pierson-Moskowitz  equation  (from  Pierson  and  Moskowitz,  1964) 


E(f)  =  ag2w5e,M”°/wfl 


provides  wave  energy  at  each  frequency  from  the  following  equation: 
where : 


w  =  2nf 


in  which/is  the  wave  frequency  in  Hertz,  g  is  gravity,  and  t/is  the  wind  speed  in  meters  per  second 
measured  at  1 9.5  m  above  the  sea  surface.  The  spectrum  E(f)  is  a  vector  of  spectral  densities  and  it  is 


a  =  0.0081 


b  =  0.74 


assumed  that  each  density  is 
the  frequency  bin. 

integrated  from  the  lower  limit  of  the  frequency  bin  to  the  upper  limit  of 

Input  Variables: 

dir 

Real 

Wave  Direction 

freql  (fieqNum) 

Real 

Beginning  Frequency  Bin  Value 

freq2  (fteqNum) 

Real 

Ending  Frequency  Bin  Value 

hsig 

Real 

Significant  Wave  Height 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 
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ifreq 

Integer 

Number  of  Frequencies  in  Input  Spectrum 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

xfrom  (dirNum) 

Real 

Direction  Array,  Direction  Wave  Energy 
Comes  From 

Output  Variables: 

esowm  (dirNum,freqNum) 

Real 

Directional  Wave  Spectrum 

Local  Variables: 

ang 

Real 

Temporary  Wave  Angle 

b 

Real 

Constant  =  0.74 

const 

Real 

Variable  in  Pierson-Moskowitz  Equation 

e 

Real 

Variable  in  Pierson-Moskowitz  Equation 

enew 

Real 

Variable  in  Pierson-Moskowitz  Equation 

gu 

Real 

Variable  in  Pierson-Moskowitz  Equation 

hs 

Real 

Set  to  Significant  Wave  Height 

hsl 

Real 

Set  to  Significant  Wave  Height 

idir 

Integer 

Direction  Loop  Counter 

ifrq 

Integer 

Frequency  Loop  Counter 

ipm 

Integer 

Set  to  1 

ratio 

Real 

Set  to  1.0 

sprd 

Real 

Directional  Spreading  Factor 

suml 

Real 

Temporary  Wave  Energy  Variable 

sum2 

Real 

Temporary  Wave  Energy  Variable 

temp 

Real 

Variable  in  Pierson-Moskowitz  Equation 

theta 

Real 

Wave  Angle 

vail 

Real 

Variable  in  Pierson-Moskowitz  Equation 

val2 

Real 

Variable  in  Pierson-Moskowitz  Equation 

wl 

Real 

Wave  Frequency  at  Beginning  of  Bin 

w2 

Real 

Wave  Frequency  at  End  of  Bin 

Subroutines  Called  from  SEAFIT  (  ): 

None. 
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SEAFIT  (  )  Called  from  Subroutines: 
WAVEFTT 

Figure  52.  Subroutine  SEAFIT  Flowchart 


START 
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Calculate  Wave 
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5.52  Subroutine  SETUP 


Subroutine  Call: 

SETUP  (pkfireq,  dl,  d2,  hrmsl,  hrms2,  etal,  kinitl,  eta2) 

Summary: 

Subroutine  SETUP  calculates  the  change  in  the  nearshore  mean  water  level  caused  by  the 
onshore  flux  of  momentum  or  the  shore-directed  Radiation  Stress.  The  presence  of  waves  causes  a 
change  in  the  total  water  depth,  which  is  defined  by  the  still  water  level  plus  the  wave-induced  set-up. 

Input  Variables: 


dl 

Real 

Corresponding  Depth 

d2 

Real 

Next  Corresponding  Depth 

etal 

Real 

Wave  Induced  Setup  at  Present  Location 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

hrms2 

Real 

Wave  Height  at  next  Onshore  Grid  Location 

kinitl 

Real 

Wave  Number 

pkfreq 

Output  Variables: 

Real 

Peak  Frequency  at  the  Center  of  the 
Frequency  Band 

eta2 

Local  Variables: 

Real 

Wave  Induced  Setup  at  New  Location 

avg_depth 

Real 

Averaged  Depth 

convrg 

Logical 

Set  to  False 

el 

Real 

Total  Average  Energy  for  Offshore  Wave 

e2 

Real 

Total  Average  Energy  for  Wave  Shoaled  and 
Refracted  Toward  the  Shore 

enl 

Real 

linear  Wave  Theory  Ratio  of  Group 

Velocity  to  Wave  Celerity 

en2 

Real 

Linear  Wave  Theory  Ratio  of  Group 

Velocity  to  Wave  Celerity 

eta_new 

Real 

Wave  Induced  Setup  Estimated  at 

New  Location 

i 

Integer 

Counter 

kl 

Real 

First  Wave  Number  Estimate 

k2 

Real 

Second  Wave  Number  Estimate 

percent_diff 

Real 

Convergence  Check 

sxxl 

Real 

Cross-Shore  Directed  Radiation  Stress 

sxx2 

Real 

Cross-Shore  Directed  Radiation  Stress 

tol 

Real 

Convergence  Tolerance 

Subroutines  Called  from  SETUP  (  ): 
WAVENUM 

SETUP  ( )  Called  from  Subroutines: 
MAIN_WAV 
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Figure  53.  Subroutine  SETUP  Flowchart 
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5.53  Subroutine  SHORTOUT 


Subroutine  Call: 

SHORTOUT  (wdir,  wspd,  j,  iimax,  dxy,  xtemp,  sural,  k,  hlmax,  h2max,  per,  pet,  thetal, 
vmax,  vmin,  width,  igamma,  bl,  it,  htemp,  wid_ii,  jgamma,  alfa,  bravo,  chrlie,  echo,  foxtrt,  golfl, 
golf2,  ihtll,  ihtl2) 

Summary: 

Subroutine  SHORTOUT  defines  the  forecasting  output  variables. 


Input  Variables: 


bl  (points) 

Real 

Bottom  Slope  Array 

dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

hlmax 

Real 

Largest  Significant  Wave  Height  in  the 

Surf  Zone 

h2max 

Real 

Largest  Maximum  Wave  Height  in  the 

Surf  Zone 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

igamma 

Integer 

Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 

iimax 

Integer 

Number  of  Calculation  Locations 

j 

Integer 

Pre-tidal  Depth  or  Still  Water  Level 

k 

Integer 

Temporary  Variable  for  Significant 

Wave  Height 

pct(4) 

Real 

Percentage  Breaker  Array 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

rk  (points,  4) 

Real 

Matrix  of  Percentage  Breakers  and 

Type  of  Breakers 

suml 

Real 

Sum  of  Wave  Length  in  the  Surf  Zone 

thetal 

Real 

Wave  Angle  at  Input  Starting  Depth 

vmax 

Real 

Maximum  Positive  Longshore 

Current  Velocity 

vmin 

Real 

Maximum  Negative  Longshore 

Current  Velocity 

wdir 

Real 

Input  Wind  Direction  -  Compass  Heading 
Wind  is  Blowing  From 

wid_ii 

Integer 

Surf  Zone  Width  Array  Index 

width 

Real 

Surf  Zone  Width 

wspd 

Real 

Input  Wind  Speed 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

Output  Variables: 
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alfa 

bravo 

chrlie 

echo 

foxtrt 

golfl 

golf2 

ihtll 

ihtl2 

jgamma 


Real 

Real 

Real 

Real 

Real 

Real 

Real 

Real 

Real 

Integer 


Local  Variables: 

11 

12 

tempi 

temp2 

xlen 


Integer 

Integer 

Real 

Real 

Real 


Subroutines  Called  from  SHORTOUT( ): 
NEW_BRK 


SHORTOUT  ( )  Called  from  Subroutines: 
CALCSURF 


Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Breaker  Angle 

Longshore  Current  Speed  and  Direction 

Number  of  Surf  Lines 

Surf  Zone  Width 

Wind  Speed 

Wind  Direction 

Temporary  Variable  Set  to 

Beach  Orientation 


Temporary  Array 
Temporary  Array 

Temporary  Variable  for  Longshore  Current 

Maximum  Calculation 

Temporary  Variable  for  Longshore  Current 

Minimum  Calculation 

Average  Wave  Length  in  Surf  Zone 
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Figure  54.  Subroutine  SHORTOUT  Flowchart 
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5.54  Subroutine  SLF_STRT 
Subroutine  Call: 


SLP_STRT  (theta,  xdelt_gr,  hnns,  per,  fqz,  fqd,  xxl,  dxy,  nnn,  eg,  hrms,  xk, 
j_ii,  M),  thetaO,  surf) 

Summary: 

Subroutine  SLF_STRT  shoals  and  refracts  waves  from  the  farthest  offshore  point  to  the 
shoreward  point  where  the  percentage  of  breaking  exceeds  the  surf  zone  criteria  of  five  percent  (5%). 
If  the  five  percent  (5%)  threshold  is  not  exceeded,  execution  halts. 


Input  Variables: 

b 

Real 

Cg 

Real 

dxy  (points) 

Real 

fqd 

Real 

fqz 

Real 

hrms 

Real 

nnn 

Integer 

per 

Real 

self_st 

Char*l 

theta 

Real 

xdelt_gr 

Real 

xk 

Real 

Output  Variables: 

Cg 

Real 

hrms 

Real 

jJi 

Integer 

10 

Real 

surf 

Logical 

thetaO 

Real 

xk 

Real 

Empirical  Factor  in  Wave  Breaking  Model 

Wave  Group  Velocity 

Corresponding  Depths  with  Tide 

Peak  Frequency  at  the  Center  of  the 

Frequency  Band 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Number  of  Points  in  Input  Depth  Array 

Peak  Period  of  Directional  Wave  Spectrum 

Self  Staring  Option  (Yes  or  No) 

Radiation  Stress  Angle 
Self-Adjusting  Cross-Shore  Grid  Step 
Wave  Number 


Wave  Group  Velocity 

Root  Mean  Square  Wave  Height 

Index  where  Wave  Probabilities 

Exceed  Threshold 

Wave  Length  Offshore  Location 

Index  Where  Percentage  of  Breakers  Is 

Exceeded  -  Start  of  Surf  Zone 

Wave  Angle  at  Grid  Offshore  Location 

Wave  Number 
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Local  Variables: 


beta 

Real 

cg2 

Real 

convg 

Real 

dp 

Real 

eb 

Real 

hrms2 

Real 

ii 

Integer 

1 

Real 

P  (4) 

Real 

rhs 

Real 

roller 

Logical 

rstart 

Real 

xkO 

Real 

Bottom  Slope 
Group  Velocity 

Convergence  Flag  (True  or  False) 

Offshore  Water  Depth 

Dissipation  Term 

Root  Mean  Square  Wave  Height 

Array  Index 

Wave  Length 

Breaker  Percentage  Array 

Right  Hand  Side  of  Energy  Equation 

Roller  Option  Flag  (True  or  False) 

Percent  Breaking  Wave  Criteria 

Offshore  Wave  Number 


Subroutines  Called  from  SLF_STRT  (  ): 

BALANCEQ 

GETJRHS 

PERCENT 

PT2 

SLF_STRT  (  )  Called  from  Subroutines: 
MAIN_WAV 
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Figure  55.  Subroutine  SLF_STRT  Flowchart 
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5.55  Subroutine  SRFSETUP 


Subroutine  Call: 

SRFSETUP  (file Jn,  file_out,  fracname,  lndname,  depname,  iyear,  imonth,  iday,  ihour,  imin, 
gamma?,  ydepth,  slope,  ydetail,  xdelt,  yrefrac,  ystr,  self_st,  hsea,  psea,  dsea,  hswell,  pswell,  dswell, 
wspd,  wdir,  tide,  spefile,  file_dat,  file_tmp,  spedepth,  file_spc,  spe_type) 

Sununary: 

Subroutine  SRFSETUP  opens  input  and  output  files.  Input  variables  are  initialized  using  data 
from  user-constructed  input  file.  The  format  of  the  input  file  is  outlined  in  Section  6.0. 

Input  Variables:  None. 


Output  Variables: 


depname 

Char*40 

Depth  Profile  File  Name 

dsea 

Real 

Input  Direction  for  Sea  Contribution 

dstart 

Real 

Input  Starting  Depth 

dswell 

Real 

Input  Swell  Direction  for  Internally 

Generated  Spectrum 

file_in 

Char*40 

Input  File  Name 

file_out 

Char*40 

Output  File  Name 

file_dat 

Char*40 

Output  Hie  Name 

file_spc 

Char*40 

Shallow  Water  Wave  Spectrum  File  Name 

file_tmp 

Char*40 

Output  Hie  Name 

fracname 

Char*40 

Wave  Refraction  File  Name 

gamma2 

Real 

Beach  Orientation,  Compass  Heading 

Directly  Toward  Beach 

gt_frg 

Integer 

Spectrum  Type 

hsea 

Real 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Pierson  Moskowitz  Spectrum 

hswell 

Real 

Input  Significant  Wave  Height  for  Internally 
Generated  Spectrum 

iday 

Integer 

Input  Day 

ihour 

Integer 

Input  Hour 

imin 

Integer 

Input  Minute 

imonth 

Integer 

Input  Month 

iyear 

Integer 

Input  Year 

lndname 

Char*40 

Input  Landing  Zone  Name 

psea 

Real 

Input  Wave  Period  for  Sea  Contribution  to 

Internally  Generated  Spectrum 

pswell 

Real 

Input  Swell  Period  for  Internally 

Generated  Spectrum 
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selfjst 

slope 

spedepth 

spefile 

spe_type 

tide 

wdir 

wspd 

xdelt 

ydepth 

ydetail 

yrefrac 

ystr 

Local  Variables: 

duml 

fend 

file_dat 

i 

iopen 

j 


Char*l 

Real 

Real 

Char*40 

Integer 

Real 

Real 

Real 

Real 

Char*l 

Char*l 

Char*l 

Char*l 


Char*80 

Integer 

Char*20 

Integer 

Integer 

Integer 


Subroutines  Called  from  SRFSETUP  ( ): 
ABORT 


SRFSETUP  (  )  Called  from  Subroutines: 
SURF 


Self  Start  Flag  (Y es  or  No) 

Bottom  Slope 

Depth  at  Offshore  Wave  Spectrum 
Selected  Wave  Spectrum  File  Name 
0  if  spefile  is  not  blank,  1  for  blank 
Input  Tide  Level 

Input  Wind  Direction,  Compass  Heading 

Wind  Blows  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Input  Depth  Profile  Used?  (Yes  or  No) 

Detailed  Output?  (Yes  or  No) 

Is  Refraction  Considered  in  Analysis? 

(Yes  or  No) 

Is  Straight  Coast  Refraction  Used?  (Yes  or  No) 


Title  Line 

File  Name  Prefix  Used  for  Building 
File  Names 

Additional  Output  File  Name 
Loop  Counter 
I/O  Status  Number 
Loop  Counter 
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Figure  56.  Subroutine  SRFSETUP  Flowchart 
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5.56  Subroutine  STRFRAC 


Subroutine  Call: 

STRFRAC  (dstart,  ifreq,  freq,  igamma,  idirec,  xfrom,  spedepth,  xcoeff,  xtheta) 

Summary: 

Subroutine  STRFRAC  calculates  wave  angle  refraction  coefficients  and  combined  shoaling 
and  refraction  coefficients  to  propagate  wave  energy  into  shallow  water. 


Input  Variables: 

dstart 

Real 

freq  (freqNum) 

Real 

idirec 

Integer 

ifreq 

Integer 

igamma 

Integer 

xfrom  (dirNum) 

Real 

Output  Variables: 

xcoeff  (dirNum, freqNum) 

Real 

xtheta  (dirNum, freqNum) 

Real 

Local  Variables: 

argl 

Real 

direc 

Real 

frd 

Real 

idir 

Integer 

ifrq 

Integer 

m 

Integer 

noprint 

Real 

shoal 

Real 

shoal2 

Real 

thetad 

Real 

thetas2 

Real 

xkd 

Real 

xk2 

Real 

Input  Starting  Depth 

Input  Wave  Spectrum  Center  Frequencies 
Number  of  Direction  Bins  in  Input  Spectrum 
Number  of  Frequency  Bins  in 
Input  Spectrum 

Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 
Direction  Array 


Wave  Height  Refraction  Coefficients 
Wave  Angle  Refraction  Coefficients 


Shallow  Water  Angle  (1)  -  Temporary 
Temporary  Direction  Angle 
Wave  Frequency 
Direction  Loop  Counter 
Frequency  Loop  Counter 
Temporary  Wave  Angle 
Wave  Component  Direction 
Temporary  Shoaling  Coefficient 
Temporary  Shoaling  Coefficient  at  Input 
Starting  Depth 

Temporary  Wave  Angle  Variable 
Temporary  Wave  Angle  Variable 
Temporary  Wave  Number  Variable 
Temporary  Wave  Number  Variable 
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xks2 


Real 


Temporary  Wave  Number  at  Input 
Starting  Depth 

xksd2  Real  Wave  Number  at  Input  Starting  Depth 

Subroutines  Called  from  STRFRAC  (  ): 

WAVENUM 

STRFRAC  ( )  Called  from  Subroutines: 

SURF 
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Figure  57.  Subroutine  STRFRAC  Flowchart 


(  START  ] 

_ '  r _ 

/  Number  of  \ 
Frequency 
Bands  in  Input 
Spectrum  Loon 


Call  WAVENUM 
at  Starting  Depth 

1 

r 

Calculate  Shoaling 
Coefficient  for  Energy 

J 

/  Number  of  \ 
Direction  Bins 
in  Input 
Spectrum  Loop 


i  i 
i  l 

Initialize  Refraction 

i  i 
t  i 

Arrays 

ii  I 

t  1  _ I _ _ _ 

I  ! 

Convert  Deep  Water 

1  ! 

Angle  to  Perpendicular  to 

I  I 

i  f 

Coast 

5.57  Subroutine  SUMMARY 


Subroutine  Call: 

SUMMARY  (spedepth,  dstart,  tide,  wspd,  wdir,  xdelt,  yrefrac,  ystr,  depname,  file_out, 
fracname,  lndname,  ydepth,  ydetail,  gamma2,  slope,  hsea,  psea,  dsea,  hswell,  pswell,  dswell,  spectra, 
spefile,  spe_type,  file_tmp) 

Summary: 


Subroutine  SUMMARY  summarizes 


documentation  and  forecaster  verification. 

Input  Variables: 

depname 

Char*40 

dsea 

Real 

dstart 

Real 

dswell 

Real 

file_out 

Char*40 

file_tmp 

Char*40 

fracname 

Char*40 

gamma2 

Real 

hsea 

Real 

hswell 

Real 

lndname 

Char*40 

psea 

Real 

pswell 

Real 

slope 

Real 

spectra 

Logical 

spedepth 

Real 

spefile 

Char*40 

spe_type 

Integer 

tide 

Real 

wdir 

Real 

wspd 

Real 

xdelt 

Real 

ydepth 

Char*l 

ydetail 

Char*l 

yrefrac 

Char*! 

the  input  information  read  to  the  output  file  for 


Depth  Profile  File  Name 

Input  Direction  for  Sea  Contribution 

Input  Starting  Depth 

Input  Swell  Direction  for  Internally 

Generated  Spectrum 

Output  File  Name  *.out 

Tempary  Output  File  Name  *.tmp 

Wave  Refraction  File  Name 

Beach  Orientation,  Compass  Heading  Directly 

Toward  Beach 

Input  Significant  Wave  Height  for  Sea 

Contribution  to  Internally  Generated  Spectrum 

Input  Significant  Wave  Height  to  Internally 

Generated  Spectrum 

Input  Landing  Zone  Name 

Input  Wave  Period  for  Sea  Contribution  to 

Internally  Generated  Spectrum 

Input  Swell  Period  for  Internally 

Generated  Spectrum 

Bottom  Slope  for  a  Constructed  Depth  Profile 

Does  Input  Spectrum  Exist?  (True  or  False) 

Wave  Input  Depth 

Selected  Wave  Spectrum  File  Name 

0  if  spefile  is  not  blank,  1  for  blank 

Input  Tide  Level 

Input  Wind  Direction  Compass  Heading  Wind 

Blows  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Input  Depth  Profile  Used?  (Yes  or  No) 

Detailed  Output?  (Yes  or  No) 

Is  Refraction  Considered  in  Analysis? 
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ystr 


Char*l 


(Yes  or  No) 

Is  Straight  Coast  Refraction  Used?  (Yes  or  No) 


Output  Variables:  None. 

Local  Variables: 

sediment  Char*40 

Subroutines  Called  from  SUMMARY  (  ): 
SUMMARY  (  )  Called  from  Subroutines: 
SURF 


Sediment  Type 
None. 
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Figure  58.  Subroutine  SUMMARY  Flowchart 
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5.58  Subroutine  SURFCAST 
Subroutine  Call: 


SURPCAST  (pet,  depname,  lndname,  slope,  ydepth,  alfa,  bravo,  chrlie,  echo,  foxtrt  eolfl 
golf2,  lhtll,  iht!2)  &  ’ 

Summary: 

Subroutine  SURFCAST  reads  input  variables  and  provides  a  short  format  summary  of  Navy 
specified  parameters.  The  subroutine  also  examines  longshore  current  direction  and  selects  the 


dominant  breaker  type. 
Input  Variables: 


alfa 

Real 

bravo 

Real 

chrlie 

Real 

depname 

Char*40 

echo 

Real 

foxtrt 

Real 

golfl 

Real 

golf2 

Real 

ihtll 

Real 

ihtl2 

Real 

lndname 

Char*40 

pet (4) 

Real 

slope 

Real 

ydepth 

Char*l 

Output  Variables: 

None. 

Local  Variables: 

foxtmp 

Real 

i 

Integer 

jdelt 

Integer 

jp(4) 

Integer 

Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Depth  Profile  File  Name 
Breaker  Angle 

Longshore  Current  Speed  and  Direction 
Number  of  Surf  Lines 
Surf  Zone  Width 

Wind  Speed  Coded  Surf  Forecast  Value 

Wind  Direction 

Input  Landing  Zone  Name 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Bottom  Slope 

Input  Depth  Profile  Used?  (Yes  or  No) 


Longshore  Current  Where  the  Sign  Indicates 

the  Direction 

Loop  Counter  Variable 

Difference  If  Any  Between  100%  and 

Sum  of  jp  (4) 

Temporary  Variable  Same  as  pct(4)  Array 


jsum 

maxp 

xmax 


Integer 

Integer 

Real 


Subroutines  Called  from  SURFCAST  (  ): 
SURFCAST  (  )  Called  from  Subroutines: 
SURF 


Check  for  Percentages  Adding  to  100% 
Indicates  Dominant  Breaker  Type 
Temporary  Variable  Used  in  Dominant  Breaker 
Type  Examination 

None. 
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Figure  59.  Subroutine  SURFCAST  Flowchart 


Output  Breaker 
Types 


Output  Littoral 
Current  Direction 


Littoral  / 
irection  / 


Output  Wind  Speed 
&  Direction  > 


5.59  Subroutine  SWLF1T 


Subroutine  Call: 

SWLFTT  (hsig,  per,  dir,  dangle,  ifreq,  idirec,  period,  esowm) 

Sununary: 

Subroutine  SWLFTT  superimposes  remotely  generated  swell  wave  energy  onto  the  existing 
directional  wave  spectrum.  The  existing  wave  spectrum  may  be  zero  or  it  may  contain  locally 
generated  sea  waves  already  added  by  the  subroutine  SEAFTT. 


Input  Variables: 


dangle 

Real 

Angle  between  Directional  Bins 

dir 

Real 

Input  Swell  Direction  for  Internally 

hsig 

Real 

Generated  Spectrum 

Significant  Wave  Height 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequency  Bins  in 

per 

Real 

Input  Spectrum 

Peak  Period  of  Directional  Wave  Spectrum 

period  (freqNum) 

Real 

Period  Array  (1  /  Frequency) 

Output  Variables: 

esowm  (dirNum, freqNum) 

Real 

Directional  Wave  Spectrum 

Local  Variables: 

dl 

Real 

Temporary  Variable  for  Distributing 
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Real 

Wave  Energy 

Temporary  Variable  for  Distributing 
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Real 

Wave  Energy 

Temporary  Variable  for  Distributing 

delt 

Real 

Wave  Energy 

Temporary  Variable  for  Distributing 

diff 

Real 

Wave  Energy 

Difference  between  Maximum  Wave 

dmin 

Real 

Period  and  Array  Value  of  Wave  Period 

Set  to  1000.0 

energy 

Real 

Swell  Energy 

ifrq 

Integer 

Frequency  Loop  Counter 

jdir 

Integer 

Swell  Direction 

jdirl 

Integer 

Direction  Bin  Index  Number 

jdir3 

Integer 

Direction  Bin  Index  Number 
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jfreq 

xdir 


Integer 

Real 


Directional  Wave  Spectrum  Wave  Number 
Wave  Direction 


Subroutines  Called  from  SWLFIT  (  ):  None. 

SWLFIT  (  )  Called  from  Subroutines: 

WAVEFTT 

Figure  60.  Subroutine  SWLFIT  Flowchart 
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5.60  Subroutine  WAVEFIT 


Subroutine  Call: 

WAVEFIT  (ifreq,  idirec,  dangle,  hsea,  psea,  dsea,  hswell,  pswell,  dswell,  freql,  fireq2,  xfrom, 
period,  esowm,  ehsig) 

Summary: 

Subroutine  WAVEFIT  initializes  the  internally  generated  directional  wave  spectrum  to  zero 
and  calls  subroutines  SEAFTT  and  SWLF1T  to  fill  the  matrix. 


Input  Variables: 

dangle 

Real 

dsea 

Real 

dswell 

Real 

freql  (freqNum) 

Real 

fireq2  (freqNum) 

Real 

hsea 

Real 

hswell 

Real 

idirec 

Integer 

ifreq 

Integer 

period  (freqNum) 

Real 

psea 

Real 

pswell 

Real 

xfrom  (dirNum) 

Real 

Output  Variables: 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

Angle  Between  Directional  Bins 
Input  Direction  for  Sea  Contribution 
to  Internally  Generated  Wave  Spectrum 
Input  Swell  Direction  for  Internally 
Generated  Spectrum 
Beginning  Frequency  Bin  Value 
Ending  Frequency  Bin  Value 
Input  Significant  Wave  Height  for  Sea 
Contribution  to  Internally  Generated 
Wave  Spectrum 

Input  Significant  Wave  Height  to 
Internally  Generated  Spectrum 
Number  of  Direction  Bins  in  Input  Spectrum 
Number  of  Frequencies  in  Input  Spectrum 
Period  Array  (1/Frequency) 

Input  Wave  Period  for  Sea  Contribution 
Input  Swell  Period  for  Internally 
Generated  Spectrum 
Direction  Array,  Direction  Wave  Energy 
Comes  From 


Significant  Wave  Height  from 
Directional  Spectrum 
Directional  Wave  Spectrum 
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Local  Variables: 


.  Integer  Direction  Loop  Counter 

^  Integer  Frequency  Loop  Counter 

Subroutines  Called  from  WAVEFIT  (  ): 

GT_SIG_H 

SEAFTT 

SWLFTT 

WAVEFIT  (  )  Called  from  Subroutines: 

GENSPEC 


Figure  61.  Subroutine  WAVEFIT  Flowchart 


[  START  ) 


(  RETURN  ] 


5.61  Subroutine  WAVENUM 


Subroutine  Call: 

WAVENUM  (fq,  dp,  xk) 

Summary: 

The  wave  dispersion  equation  is  solved  for  the  wave  number  through  numerical  iteration.  A 
relative  change  of  less  than  .0005  is  required  and  the  maximum  number  of  iterations  is  150.  If 
convergence  is  not  obtained  within  150  iterations,  a  shallow  water  approximation  is  employed. 

Input  Variables: 

dp  Real  Offshore  Water  Depth 

fq  Real  Wave  Frequency 

Output  Variables: 

xk  Real  Wave  Number 

Local  Variables: 

const  Real  Shallow  Water  Criteria  Constant 

diff  Real  Percent  Difference  between  Wave 

Number  Estimates 

est  Real  Estimate  of  Wave  Number 

I  Integer  Loop  Counter 

it  Integer  Loop  Limit  -  Set  to  150 

Subroutines  Called  from  WAVENUM  (  ):  None. 

WAVENUM  ( )  Called  from  Subroutines: 

INITLIZE 

PT2 

RADJST1 

RAD_ST2 

SETUP 

STRFRAC 
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Figure  62.  Subroutine  WAVENUM  Flowchart 


[  START 


_ I 

Set  Maximum  Number 
of  Iterations 
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5.62  Subroutine  WEIGHTFN 


Subroutine  Call: 


WEIGHTFN  (dp,  hrms,  h,  w_h) 


Summary: 

Subroutine  WEIGHTFN  calculates  the  weighting  function  used  to  describe  the  distribution  of 


breaking  waves  across  the  surf  zone. 


Input  Variables: 

dp 

h 

hrms 

Output  Variables: 
w_h 

Local  Variables: 
m 

temp 

tol 


Real 

Real 

Real 


Real 


Real 

Real 

Real 


Offshore  Water  Depth 
Wave  Height 

Root  Mean  Square  Wave  Height 


Output  Weighting  Function 


Multiplier 
Weighting  Function 
Set  to  -700.0 


Subroutines  Called  from  WEIGHTFN  (  ): 
GET_M 

WEIGHTFN  (  )  Called  from  Subroutines: 
F2 
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Figure  63.  Subroutine  WEIGHTFN  Flowchart 


[  START  ) 

J 


1 


Calculate 

Weighting 

Function 


[ 

[  RETURN  "[ 
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5.63  Subroutine  ZONE1 


Subroutine  Call: 

ZONE1  (jji,  iimax,  dxy,  xtemp,  htemp,  ptemp,  thetatemp,  xktemp,  v,  distmax,  vmax,  vmin, 
thetamin,  thetamax,  suml,  width,  j,  k,  hlmax,  h2max,  wid_ii) 

Summary: 

Subroutine  ZONE1  calculates  the  preliminary  surf  forecast  values  and  surf  zone  parameters. 


Input  Variables: 

distmax 

Real 

Farthest  Distance  Offshore 

dxy  (points) 

Real 

Pre-Tidal  Depth  or  Still  Water  Level 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

iimax 

Integer 

Number  of  Calculation  Locations 

j-ii 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

ptemp  (points) 

Real 

Percentage  of  Breaking  Waves  and 

Breaker  Types 

v  (points) 

Real 

Longshore  Current 

xktemp  (points) 

Real 

Temporary  Variable  for  Wave  Number 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

Output  Variables: 

hlmax 

Real 

Maximum  Significant  Wave  Height 

h2max 

Real 

Maximum  Wave  Height 

j 

Integer 

Array  Index  Where  Maximum  Significant 
Wave  Height  Occurs 

k 

Integer 

Temporary  Variable  Number  of  Points  in 
Cross-Shore  Transect 

suml 

Real 

Summation  of  Wave  Lengths  Across  the 

Surf  Zone 

vmax 

Real 

Maximum  Positive  Longshore  Current 
Velocity 

vmin 

Real 

Maximum  Negative  Longshore 

Current  Velocity 

wid_ii 

Integer 

Array  Index  for  X- value  at  Surf 

Zone  Boundary 

width 

Real 

Surf  Zone  Width 
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Local  Variables: 


dpi 

hdep 

hmax 

houtl 

hrmsl 

ii 

vingl 

wlen 

xoffl 


Real 

Real 

Real 

Real 

Real 

Integer 

Real 

Real 

Real 


Subroutines  Called  from  ZONE1  (  ): 
ZONE1  (  )  Called  from  Subroutines: 
CALCSURF 


Offshore  Depth  in  Feet 
Limiting  Breaking  Depth 
Temporaiy  Variable  for 
Maximum  Wave  Height 
Temporary  Variable  for  Significant 
Wave  Height 

Root  Mean  Square  Wave  Height 
Loop  Index 

Longshore  Current  Velocity  in  Knots 
Wave  Length 
Distance  Offshore 

None. 
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Figure  64.  Subroutine  ZONE1  Flowchart 
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5.64  Function  F2 


Function  Call: 

F2  (h,  hrms,  dp,  p_flag) 

Summary: 

Function  F2  evaluates  the  Rayleigh  probability  distribution  function  for  a  given  wave  height 
value,  for  a  selected  weighting  function. 


Input  Variables: 

dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

hrms 

Real 

Root  Mean  Square  Wave  Height 

P_flag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

Output  Variables: 

f2 

Real 

Weighted  Rayleigh  Distribution 

Local  Variables: 

P-h 

Real 

Rayleigh  Probability  Distribution 

temp 

Real 

Exponent  Term  in  Rayleigh  Distribution 

tol 

Real 

Tolerance  Value  Set  to  -700.0 

w_h 

Real 

Weighting  Function 

Subroutines  Called  from  F2  (  ): 
WEIGHTFN 

F2  ( )  Called  from  Functions: 


INTEGRAT 


Figure  65.  Function  F2  Flowchart 


Calculate  Weighted 
Rayleigh  Distribution 
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5.65  Function  F3 


Function  Call: 


F3  (farms,  theta,  Cg,  dp,  meanjfreq,  xk,  roller) 

Summary: 


Function  F3  returns  values  for  the  LHS  of  the  energy  equation. 
Input  Variables: 


Cg 

dp 

hrms 

meanjfreq 

roller 

theta 

xk 

Output  Variables: 
f3 

Local  Variables: 


Real 

Real 

Real 

Real 

Logical 

Real 

Real 


Real 


e_roller  Real 

e_wave  Real 

Subroutines  Called  from  F3  ( ): 

CALCROLL 

GET_WAVE 


Wave  Group  Velocity 
Offshore  Water  Depth 
Root  Mean  Square  Wave  Height 
Directional  Spectrum  Value 
Roller  Option  Flag  (True  or  False) 
Wave  Angle 
Wave  Number 


Total  Energy 


Roller  Contribution  to  the  Energy  Equation 
Wave  Contribution  to  the  Energy  Equation 


F3  (  )  Called  from  Subroutines: 
BALANCEQ 


Figure  66.  Function  F3  Flowchart 
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5.66  Function  INTEGRAT 


Function  Call: 

INTEGRAT  (xo,  xn,  hrms,  dp,  p_flag) 
Summary: 


Function  INTEGRAT  evaluates  an  integral  numerically  using  the  trapezoidal  rule.  Function 
{F2}  is  called  to  evaluate  the  integral  at  upper  and  lower  limits.  The  function  applies  the  trapezoidal 
integration  method  to  estimate  the  wave  height  at  a  particular  depth  from  a  weighted  distribution. 
Input  Variables: 


dp 

hrms 

P_flag 

xn 

xo 

Output  Variables: 
integrat 


Real 

Real 

Logical 

Real 

Real 


Real 


Local  Variables: 

delt 

f_xn 

f_xo 

f2 


Real 

Real 

Real 

Real 


1  Integer 

numit  Integer 

sum  Real 

Real 

Functions  Called  from  INTEGRAT  ( ): 


F2 


Farthest  Offshore  Water  Depth 
Root  Mean  Square  Wave  Height 
Weighting  Factor  Flag  (True  or  False) 
Upper  Limit  of  Integration  =  5  *  hrms 
Lower  Limit  of  Integration  =  0.0 


Wave  Height  Distribution  Calculated  for  a 
Specific  Location 


Step  Between  Intervals 
f(x)  Evaluated  at  Upper  Limit 
f(x)  Evaluated  at  Lower  Limit 
Wave  Height  Distribution 
Weighting  Function 
Loop  Variable 

Set  to  100  -  Number  of  Iterations  Examined 
Over  Integral 

Summary  Results  from  Function  F2 
Integration  Step  Location 
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INTEGRAT  (  )  Called  from  Subroutines: 

CALCHB3 

PERCENT 

Figure  67.  Function  INTEGRAT  Flowchart 
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5.67  Include  File:  COMMON.INC 
Summary: 


The  include  file  COMMON.INC  contains  all  the  parameters  set  for  the  SURF  Model.g 


Defined  Parameters: 

deal 

Real 

0.3048  -  Feet  to  Meters  Conversion 

degrad 

Real 

PI  /  180.0  -  Conversion  from 

Degrees  to  Radians 

dirNum 

Integer 

180  -  Array  Dimension  Used  for 

Direction  Arrays 

freqNum 

Integer 

50  -  Array  Dimension  Used  for 

Frequency  Arrays 

g 

Real 

9.8 

gamma 

•  •  * 

Real 

0.42  -  Empirical  Wave  Height  Factor 

lumt 

Integer 

Output  File  Unit 

pi 

Real 

3.14159265 

points 

Integer 

500  -  Array  Dimension  Used  for  all  Input 

raddeg 

Depth  Arrays 

Real 

180.0  /  pi  -  Conversion  from 

Radians  to  Degrees 

riio 

Integer 

1030  -  Water  Density 

riioair 

Real 

1.2  -  Air  Density 

sigma 

Real 

sigma_deg  *  degrad 

sigma-deg 

Real 

5.0  -  Angle  in  Degrees  between  Wave  and 

tpi 

Real 

Roller  in  the  Thomton/Lippman  Model  (1996) 

2*3.14159265 

zone_pct 

Real 

10%  Surf  Zone  Width  Percent  of  Breaking 
Waves 
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APPENDICES 


Appendix  A.  INPUT  AND  OUTPUT  DATA  FORMATS 
1.  Input  File  Formats  and  model  options 

This  section  gives  the  formats  for  the  files  read  or  produced  by  SURF. 
1.1  SURF  Input  File 


The  SURF  input  file  contains  12  lines.  Some  of  the  lines  may  be  blank;  some  are  required.  The 
format  for  each  line  of  the  input  file  is  al  follows: 


Line  Description 

1  Input  File  Name 


Type  Units  Range 

Char* 40  -  - 


The  entry  in  line  1  must  be  the  exact  name  of  the  input  file.  The  first  character  of  the  file  name  must  be 
in  column  1.  The  file  name  is  limited  to  40  characters. 


Line  Description 

2  Date  and  Time  YYYYMMDDHH 


Type  Units  Range 

Char*10  -  - 


Line  2  is  date-time  information  in  the  form  YYYYMMDDHH.  SURF  simply  reads  this  line  and  prints 
it  out  in  the  output  file. 


Line  Description 

3  Landing  Zone  Name 


Type  Units  Range 

Char*40  -  - 


Line  3  is  a  description  of  the  beach.  The  string  in  line  3  cannot  be  longer  than  40  characters  or  the 
string  will  be  truncated.  This  line  can  be  blank,  but  no  information  to  identify  the  beach  will  appear  in 
the  output  file. 


Line 

4 


Description  Type 

Input  Depth  Profile  File  Name  Char*40 


Units 


Range 


Line  4  is  the  name  of  the  input  depth  profile.  The  depth  profile  file  name  is  limited  to  40  characters. 


Description  Type  Units  Range 

Sediment  Type  Integer  -  1-10 


An  entry  in  line  5  must  be  given  if  no  depth  profile  file  is  included  in  line  4.  If  a  depth  profile  is 
specified  in  line  4,  this  line  should  be  left  as  blank.  Allowable  entries  for  bottom  composition  are  as 
follows 

1  =  Boulders 

2  =  Cobble 

3  =  Pebbles 

4  =  Granules 

5  =  Very  Coarse  Sand 

6  =  Coarse  Sand 

7  =  Medium  Sand 
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Line 

6 


8  =  Fine  Sand 

9  =  Very  Fine  Sand 

10  =  Silt 


Units 

Degrees 


Range 

0-359 


Description 

Compass  Heading  Towards  Beach 


Type 

Real 


The  compass  heading  toward  the  beach  is  the  direction  from  sea  to  beach,  perpendicular  to  the  beach. 
Some  examples  of  beach  orientation  are  shown  in  Fig.  1,  part  (a). 


Line  Description 

7  Wave  Input  Depth 


Type  Units  Range 

Real  Feet  >  0 


Line  7  is  the  depth  in  feet  at  the  location  of  the  input  waves.  The  input  waves  can  be  in  two  formats: 

(1)  a  directional  wave  spectrum  from  a  file  given  in  line  8.  Straight  coast  refraction  will  be  applied 
if  the  depth  is  deeper  than  available  depth  profile.  If  line  9  (wave  refraction  file)  is  not  blank, 
this  depth  corresponds  to  the  output  depth  where  transformation  coefficients  are  applied  to 
offshore  input  wave.  Further  illustration  is  included  in  the  section  1.5. 

(2)  sea  and  swell  parameters  in  line  10,  which  are  used  to  generate  a  synthetic  directional  wave 
spectrum  within  SURF; 


Line  Description  Type  Units  Range 

8  Input  Wave  Spectrum  File  Name  Char*40  -  - 

Line  8  is  the  name  of  the  optional  input  directional  wave  spectrum  file.  If  a  file  is  entered  here  then  any 
wave  input  information  line  10  is  ignored  during  SURF  execution. 

Line  Description  Type  Units  Range 

9  Input  Wave  Refraction  File  Name  Char*40  -  - 

Line  9  is  the  name  of  the  input  refraction  and  shoaling  file.  It  should  be  noted  that  the  depth  at  the 
offshore  boundary  of  the  wave  refraction  computation  domain  should  be  the  same  as  offshore  wave 
spectrum  input  depth.  A  wave  spectrum  from  line  8  or  wave  input  from  line  10  will  be  modified  by 
the  refraction  angles  and  shoaling  coefficients  in  this  file.  If  this  line  is  blank,  then  simple  refraction 
and  shoaling  based  on  a  straight  coast  assumption,  i.e.  parallel  bottom  contours,  will  be  applied. 


Line 

Description 

Type 

Units 

Range 

10 

Sea  Wave  Height 

Real 

Feet 

>  0 

Sea  Wave  Period 

Real 

Seconds 

.  1  -  30 

Sea  Wave  Direction 

Real 

Degrees 

0  -  359 

Swell  Wave  Height 

Real 

Feet 

>  0 

Swell  Wave  Period 

Real 

Seconds 

1-30 

Swell  Wave  Direction 

Real 

Degrees 

0  -  359 

Wave  direction  is  the  direction  from  which  waves  come  in  degrees  from  North.  Some  examples  of 
wave  direction  are  shown  in  Fig.  1 ,  part  (b).  If  no  directional  wave  spectrum  file  is  given  in  line  8  then 
the  model  will  produce  a  directional  wave  spectrum  based  on  the  sea  and  swell  parameters  given  in 
this  line.  If  a  refraction-shoaling  file  is  included  then  the  internally  generated  spectrum  will  be 
refracted  and  shoaled  to  the  depth  in  line  7. 


Type 

Real 

Real 


Line 

11 


Description 
Wind  Speed 
Wind  Direction 
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Units 

Knots 

Degrees 


Range 

>  0 

0  -  359 


Tide  Elevation 


Real 


Feet 


+  or 


Une  1  i  gives  wind  and  tide  information.  Wind  direction  is  the  direction  from  which  wind  comes 
in  degrees  from  North. 


Line  Description 

12  Output  Grid  Spacing 


Type  Units  Range 

Real  Feet  see  note 


An  entry  must  be  made  in  line  12.  If  line  12  is  negative  then  a  short  output  will  be  produced. 

Note:  the  range  of  intervals  is  limited  by  array  sizes  and  by  the  surf  zone  width  computed  by  the 
model.  Error  messages  will  warn  the  user  if  the  intervals  are  too  small,  say  less  than  2  ft,  or  too  large. 


270° 


180® 


77777777777777777777 

Beach 


(a)  Beach  Orientation  Definition.  Arrows  show  sight  lines  from  deep 
water  toward  Ihe  beaches. . 


North 


Wave  Direction  0®  Wave  Direction  90® 


Waves 


Wave  Direction  270® 

- - - — . 


(b)  Wave  Diration  Definition.  Directions  are  those  from  which  waves 
comein  degrees  r&thr£  to  North. 


Fig.  1  Beach  orientation  and  wave  direction  definitions. 
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1.2  Depth  File 


Type 
Char* 80 


Range 


Line 

1 


Description 

Title 


line  1  is  a  simple  alphanumeric  identifier.  The  information  in  this  line  is  not  used  in  SURF. 

Line  Description  Type  Range 

2  Units  for  Distance  Offshore  Integer  1,2  or  3 

Line  2  identifies  the  units  of  offshore  distances  associated  with  the  entries  in  line  4  and  after. 

1  -  Distances  in  Feet 

2  -  Distances  in  Meters 

3  -  Distances  in  Yards 

Line  Description  Type  Range 

3  Units  for  Depth  Integer  1,2  or  3 


line  3  identifies  the  units  of  the  depths  associated  with  the  entries  in  line  4  and  after. 

1  -  Depths  in  Feet 

2  -  Depths  in  Meters 

3  -  Depths  in  Fathoms 


Line  Description 

4+  Index  number 

Distance  offshore 
Depth  positive  down 


Type  Range 

Integer  1  -  500 

Real  - 

Real  - 


The  depth  profile  is  contained  in  lines  4  and  after.  The  distance  coordinate  is  zero  at  the  water’s  edge 
and  increases  offshore.  The  depths  associated  with  each  distance  are  positive  down.  See  Appendix  B 
for  a  sample  input  depth  profile  file. 


1.3  Directional  Wave  Spectrum  File 


The  input  directional  wave  spectrum  file  contains  nine  preliminary  lines  of  information  followed  by 
blocks  of  data,  where  each  block  is  associated  with  a  frequency  band.  The  elements  of  each  block  are 
values  of  spectral  energy  density  in  units  of  meters-squared  per  hertz  per  radian. 

Lines  1-3  identify  the  time  and  location  of  the  spectrum.  This  information  is  not  used  by  the  model  in 
calculating  wave  or  surf  parameters. 


Line 

Description 

Type 

Units 

Range 

1 

Longitude 

Real 

Degrees 

-180  -  180 

2 

Latitude 

Real 

Degrees 

-90  -  90 

3 

Date  -  (YYYYMMDD) 

Real 

— 

Line 

Description 

Type 

Units 

Range 

4 

Number  of  Angles 

Integer 

1  -  180 

Line  4  gives  the  number  of  direction  bins  in  the  directional  wave  spectrum.  The  number  in  line  4  must 
equal  the  number  of  rows  times  the  number  of  columns  in  line  5. 
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Line 

5 


Description 
Number  of  Rows 
Number  of  Columns 


Type 

Integer 

Integer 


Units 


Range 

+  number 
+  number 


TTiis  line  gives  information  for  reading  each  block  of  spectral  energy  densities.  Each  block  has  the 
same  number  of  elements,  which  is  the  number  of  rows  times  the  number  of  columns.  Note  that  the 
number  of  elements  must  be  an  even  number.  If  the  input  directional  wave  spectrum  has  24  direction 
bins  then  acceptable  pairs  of  row-column  combinations  are :  24  1;  12  2;  6  4;  3  8;  8  3;  4  6;  2 12;  1 24. 


Line  Description  Type  Units  Range 

6  Number  Frequency  Bands  Integer  -  1-50 

Line  6  contains  the  number  of  frequency  bins  in  the  directional  wave  spectrum. 

Line  Description  Type  Units  Range 

7  Initial  Direction  Real  Degrees  0  —  359 


The  directional  bands  associated  with  the  spectrum  must  increase  monotonically.  Line  7  gives  the 
initial  direction,  which  will  be  the  smallest  angular  value  in  degrees,  positive  clockwise  from  North. 

Line  Description  Type  Units  Range 

8  Width  of  Direction  Bin  Real  Degrees  2  —  180 

The  number  of  directional  bands  is  given  in  line  8. 


Note:  the  width  of  the  direction  bins  in  degrees  times  the  number  of  direction  bins  must  equal  360 
degrees. 

Line  Description  Type  Units 

9  Direction  of  Waves  Integer  - 

1  -  Direction  waves  are  coming  from 

2  -  Direction  waves  are  going  to 

Following  the  initial  nine  lines,  are  blocks  of  values  of  spectral  energy  density  in  units  of  meters- 
squared  per  hertz  per  radian.  The  first  line  of  each  block  will  contain  the  lower,  center  and  upper 
frequency  of  the  frequency  band  associated  with  that  block.  The  block  of  values  is  a  rectangular  matrix 
of  values  in  order  from  left  to  right  being  from  left  to  right  in  direction  in  increments  of  the  directional 
bandwidth  given  in  line  8.  The  block  of  data  must  represent  directions  covering  360  degrees  from  the 
initial  directional  clockwise.  In  general,  the  format  of  each  block  is  a  follows: 

Directional  Wave  Spectrum  -  Blocks  are  repeated  for  each  Frequency  Bin 
Line  Description  Type 

10  Bin  Number  Integer 

Lower  Limit  of  Frequency  Bin  Real 

Center  of  Frequency  Bin  Real 

Upper  Limit  of  Frequency  Bin  Real 

11+  Directional  Wave  Spectrum  Real 

The  elements  of  each  block  of  values  comprising  the  spectral  energy  densities  for  a  given  frequency  are 
in  the  form  of  a  rectangular  matrix  of  numbers  of  the  number  of  rows  times  the  number  of  columns,  as 
in  line  5. 


Units 

Range 

1-50 

hertz 

>  «  0 

hertz 

>  =  0 

hertz 

>  «  0 

m2/Hz/rad 

>  =  0 

Range 

1  or  2 
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1.4  Input  Wave  Refraction  and  Shoaling  Input  File 

Using  the  input  wave  refraction  and  shoaling  input  file  is  an  advanced  procedure.  The  refraction  and 
shoaling  files  used  to  modify  an  input  directional  wave  spectrum  to  a  spectrum  representative  of 
conditions  at  the  depth  given  in  line  7  of  the  SURF  input  file. 

Line  Description  Type  Units  Range 

1  Header  Character  -  - 

2  Header  Character  -  - 

3  Input  and  Output  Real  Feet  - 

Depths 

Lines  1-3  are  strings  of  identifying  text.  The  information  is  not  used  in  computation.  In  line  3,  input 
depth  is  the  offshore  boundary  depth,  and  output  depth  corresponds  to  the  depth  where  the 
transformation  coefficients  are  saved,  i.e.  the  spedepth  of  line  7  of  surf  input  file. 


Line  Description  Type  Units  Range 

4  Number  of  Angles  Integer  -  1-180 

5  Number  of  Rows  Integer  -  +  number 

Number  of  Columns  Integer  -  +  number 

6  Number  of  Freq.  Bins  Integer  -  1-50 

7  Initial  Direction  Real  Degrees  0  -  359 

8  Width  of  Direction  Bin  Real  Degrees  2  -  180 

9  Direction  of  Waves  Integer  -  1  or  2 

1  -  Direction  waves  are  coming  from 

2  -  Direction  waves  are  going  to 


Lines  4-9  are  similar  to  those  in  the  input  directional  wave  spectrum  file. 

Refraction  Angles  -  This  section  is  repeated  for  each  frequency 


Line  Description  Type  Units  Range 

10  Bin  Number  Integer  -  1-50 

Lower  Limit  of  Frequency  Bin  Real  Hertz  >  =  0 

Center  of  Frequency  Bin  Real  Hertz  >  -  0 

Upper  Limit  of  Frequency  Bin  Real  Hertz  >  =  0 


11+  Refraction  Angles  Real  Degrees  0  -  359 

End  of  Refraction  Angles 

The  elements  of  each  block  of  values  comprising  the  refraction  angles  for  a  given  frequency  are  in 
the  form  of  a  rectangular  matrix  with  the  number  of  rows  and  columns  in  line  5.  Pad  fields  with 
zeros,  if  necessary. 

Line  Description  Type  Units  Range 

Line  A+l  Header  1  for  Shoaling  Coefficients  Char*80  -  - 

Line  A+2  Header  2  for  Shoaling  Coefficients  Char*80  -  - 

Line  A+3  Header  3  for  Shoaling  Coefficients  Char*80  -  - -- 

The  Line  a+  numbering  above  and  below  denotes  information  after  the  block  of  refraction  angles. 

Shoaling  Coefficients  -  This  section  is  repeated  for  each  frequency 


Line 

Description 

Type 

Units 

Range 

Line 

A+4  Bin  Number 

Integer 

1-50 

Lower  Limit  of 

Freq 

Bin 

Real 

Hertz 

>  =  0 

Center  of  Freq 

Bin 

Real 

Hertz 

>  =  0 

Upper  Limit  of 

Freq 

Bin 

Real 

Hertz 

>  =  0 
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Line  A+5+  Shoaling  Coefficients  Real 
End  of  Shoaling  Coefficients 


m2/m2 


The  elements  of  each  block  of  values  comprising  the  shoaling  coefficients  for  a  given  frequency  are  in 
the  form  of  a  rectangular  matrix  of  values  with  the  number  of  rows  and  columns  given  in  line  5.  Pad 
fields  with  zeros,  if  necessary. 

Note:  The  angles  and  coefficients  in  this  file  must  be  defined  over  the  entire  range  (0, 360)  degrees.  A 
partial  sector  definition  (e.g.  0  to  180  degrees)  will  cause  errors.  If  the  input  data  are  not  available  over 
the  entire  range  pad  the  refraction  and  direction  bins  with  zeros. 

1.5  Model  Options 

This  section  gives  options  in  SURF  that  control  wave  refraction,  equilibrium  profile  option  and  the 
output  files. 

Wave  Refraction  Options 

In  general,  the  depth  profile  should  cover  depths  to  around  30  ft.  If  the  depth  of  input  waves  is 
deeper  than  the  deepest  depth  in  the  profile,  i.e.  outside  of  the  profile-covered  area,  two  options  are 
available  to  consider  wave  refraction  to  bring  the  input  waves  to  the  edge  of  the  area  over  which  SURF 
is  to  operate.  If  no  bathymetry  information  is  available,  straight  coast  refraction,  assuming  parallel 
bottom  contours,  will  be  used.  If  bathymetry  information  available,  one  can  use  the  wave  modeling 
option  where  transformation  coefficients  for  refraction  and  shoaling  are  computed.  This  option  is 
generally  only  used  when  the  bathymetry  is  complicated. 


As  illustrated  in  Fig.  2,  straight  coast  refraction  brings  the  wave  input  to  the  edge  of  the  profile-covered 
area  starting  point.  It  should  be  noted  that  if  spedepth  is  inside  the  profile-covered  area,  then  it 
becomes  starting  point.  Then  no  additional  wave  refraction  will  be  applied.  The  first  output  point 
corresponds  to  a  location  where  the  percent  of  wave  breaking  has  reached  5%.  This  avoids  a  long 
listing  of  surf  output  over  long  stretches  of  flat,  gently  sloping  bottoms. 
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Fig.  2  Illustration  of  straight  coast  refraction  option.  Straight  coast  refraction  brings  wave  input 
to  the  edge  of  the  profile-covered  area,  i.e.  the  starting  point. 

For  the  wave  modeling  option,  wave  models  such  as  REFDIF,  STWAVE  and  SWAN  pre-compute 
needed  transformation  coefficients  for  a  given  bathymetry.  Input  line  9  specifies  the  wave  refraction 
file.  As  illustrated  in  Fig.  3,  the  spedepth  (input  line  7)  corresponds  to  the  output  depth  of  the 
refraction  computation.  For  accuracy,  it  is  requires  that  the  output  depth  is  within  the  profile-covered 
area.  This  is  because  no  additional  straight  coast  refraction  will  be  applied  if  it  falls  outside  of  the 
profile-covered  area.  Hie  output  depth  should  not  be  too  shallow  (e.g.  within  the  surf  zone),  because 
the  transformation  coefficient  approach  assumes  that  no  depth  induced  wave  breaking  has  occurred  at 
the  output  point.  It  is  recommended  that  the  output  depth  should  be  around  25  to  30  ft  or  deeper 
depending  on  the  bathymetry  and  wave  conditions.  The  offshore  wave  input  location  needs  to  be  at  the 
same  depth  as  the  offshore  boundary  of  the  refraction  file  computation. 
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Fig.  3  Illustration  of  wave  modeling  option. 

Equilibrium  Profile  Option 


The  equilibrium  profile,  based  on  sediment  size,  is  used  if  a  depth  profile  is  not  available.  In  the 
code,  its  maximum  depth,  also  the  starting  depth,  is  set  to  10  m,  except  for  the  wave  refraction  file 
option  in  which  the  maximum  depth  corresponds  to  spedepth. 

Wave  Spectrum  Output  Option 

To  obtain  an  output  directional  wave  spectrum  file,  place  the  character  “+”  in  front  of  the 
directional  wave  spectrum  file  name  in  line  9  of  the  basic  input  file.  The  output  file  will  give  the 
directional  wave  spectrum  associated  with  the  output  point,  i.e.  the  spedepth  depth.  The  output 
spectrum  file  will  have  the  same  file  name  as  the  input  file  name  but  with  the  extension  .  dws. 

Short  Output 

In  SURF,  the  user  can  control  the  amount  of  data  in  the  output  file.  If  line  12  contains  a  zero  or  a 
negative  number,  a  short  output,  without  cross  shore  profiles  of  surf  parameters,  will  be  produced.  The 
short  output  is  similar  in  format  to  naval  surf  observations. 
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2.  Output  File  Formats 


2.1  Basic  Output  File 

The  SURF  detailed  output  has  three  output  sections  delineated  by  lines  of  asterisks.  The  first 
section  contains  input  parameters  describing  the  directional  wave  spectrum.  The  second  section  is  the 
coded  surf  forecast  with  variables  specific  to  military  surf  observations.  The  final  section  is  the 
optional  detailed  surf  output,  which  is  comprised  of  a  table  of  cross  shore  surf  zone  parameter.  These 
parameters  include  cross  shore  distance,  depth,  wave  height,  wave  breaking,  wave  angle  and  longshore 
current.  The  filename  generated  has  the  same  name  as  the  input  file  but  the  extension  is  .  out. 


Section  1 

Line 

Description 

Type 

Units 

Line  1 

Surf  Forecast  Header 

Character 

Line  2 

Blank  Line 

Line  3 

SURF  Model  Version 

Character 

Line  4 

Date  and  Time  of  Forecast 

Character 

Line  5 

Output  File  Name  Information 

Character 

Line  6 

Landing  Zone  Name 

Character 

Line  7 

Sight  Line  Toward  Beach 

Real 

Degrees 

Line  8 

Depth  Profile  Name  or  Beach  Sediment 

Type  Character 

Line  9 

Wave  Input  Depth 

Real 

Feet 

Line  10 

Spectrum  Usage  Text 

Character 

or 

Sea  Wave  Height 

Real 

Feet 

Sea  Period 

Real 

Seconds 

Sea  Direction 

Real 

Degrees 

Line  11 

Spectrum  File  Name 

Character 

or 

Swell  Wave  Height 

Real 

Feet 

Swell  Period 

Real 

Seconds 

Swell  Direction 

Real 

Degrees 

Line  12 

Wind  Speed 

Real 

Knots 

Line  13 

Wind  Direction 

Real 

Degrees 

Line  14 

Tide  Level 

Real 

Feet 

Line  15 

Blank  Line 

Character 

Line  16 

Wave  Refraction  Option 

Character 

- Line 

17  Starting  Depth 

Real 

Feet 

Line  18 

Output  Interval 

Real 

Feet 

Line  19 

Computational  grid  Spacing 

Real 

Feet 

Line  20 

Input  Spectrum  Type 

Character 

Line  21 

Significant  Wave  Height  Offshore 

Real 

Feet 

Line  22 

Wave  Peak  Period 

Real 

Seconds 

line  23 

Average  Wave  Direction 

Real 

Degrees 

line  24 

Percent  Breaking  Waves  at  Starting  Depth 

Real 

Percent 

It  should  be  noted  that  starting  depth  on  line  17  is  the  depth  after  offshore  waves  have  brought  to  the 
edge  of  the  profile-covered  area  through  either  straight  coast  refraction  or  refraction  file  computation. 
This  depth  depends  on  the  depth  profile,  tide  and  wave  input  (spedepth)  location. 

Section  2 


Line 


Description 
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Type 


Units 


Iiine 

1 

Line 

2 

Line 

3 

Line 

4 

Line 

5 

Line 

6 

Line 

7 

Line 

8 

Line 

9 

Line 

10 

Line 

11 

Line 

12 

Line 

13 

Line 

14 

Line 

15 

Section  3 

Line 

Line 

1 

Line 

2 

Line 

3 

Line 

4 

Line 

5 

Line 

6 

Line 

7 

Line 

8- 

Code  Surf  Forecast 
Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Dominant  Breaker  Type 
Breaker  Percentages 
Breaker  Angle 
Littoral  Current 
Number  of  Surf  Lines 
Surf  Zone  Width 
Wind  Speed 
Average  wave  length 
Wind  Direction 
Blank  Line 
Modified  Surf  Index 


Description 
Blank  Line 

Heading  -  Detailed  Surf  Output 
Blank  Line 
Text  Heading  Line 
Text  Heading  Line 
Text  Heading  Line  -  Units 
Blank  Line 
EOF  Index  Number 

Distance  Offshore 
Water  Depth 

Significant  Breaker  Height 
Maximum  Breaker  Height 
Percent  Breaking  Waves 
Breaker  Angle 
Littoral  Current 


Character 

Real 

Feet 

Real 

Feet 

Real 

Seconds 

Character 

Character 

Percent 

Real 

Degrees 

Real 

Knots 

Real 

Real 

Feet 

Real 

Knots 

Real 

Feet 

Real 

Degrees 

Character 

Real 

Type 

Units 

Character 

Character 

Character 

Character 

Character 

Character 

Integer 

Real 

Feet 

Real 

Feet 

Real 

Feet 

Real 

Feet 

Real 

Percent 

Real 

Degrees 

Real 

Knots 

The  first  output  point  in  line  8  in  section  3  corresponds  to  a  point  where  percent  of  wave  breaking 
has  reached  5%.  6 


2.2  Data  Only  Output  File 


The  data  only  output  file  contains  the  same  information  in  the  same  format  as  the  section  3  of  the 
detailed  model  output,  except  the  file  does  not  contain  header  information.  It  is  useful  in  graphic 
applications. 


23  Shallow  Water  Directional  Wave  Spectrum 


The  shallow  water  directional  wave  spectrum  output  file  is  created  when  the  first  character  of  line  6 
in  the  basic  input  file  is  a  This  file  has  the  same  file  name  as  the  input  file  except  that  the  file 
extension  will  be  .  dws.  The  first  row  contains  the  center  frequencies  of  the  directional  wave  energy 
spectrum.  The  first  column  defines  the  wave  directions  of  the  directional  wave  energy  spectrum.  The 
remaining  matrix  elements  comprise  the  directional  wave  energy  spectrum. 


Description 

1  Frequency  Bins 

Column  1  Wave  Direction 

Other  elements  Spectral  Energy  Density 


Type 

Real 

Real 

Real 


Units  Range 

hertz  0-0.5 

degrees  0-359 

m2/(Hz-rad)  0  -  999 


180 


Appendix  B.  Error  Message  Descriptions 


Error  Message 

Subroutine 

Generating 

Error 

Suggested  Solution  to  Resolve 
Error 

| 

Error  1 15  -  Opening  Directional  Wave 
Spectrum  File. 

readspec 

Check  wave  spectrum  file  name  in 
the  input  file- line  5.  Verify  the 
location  of  the  spectrum  file  is  the 
same  as  the  input  file. 

Error  120  -  Opening  input  file. 

srf setup 

Check  the  name  of  the  input  file 
typed  at  the  command  prompt 
(surf  32  <  f  n .  in)  or  the  name 
typed  during  execution  (Enter 
f  n .  in). 

Error  125  -  Opening  of  Input  Depth  File. 

c_in_dep 

Check  depth  profile  file  name  in 
the  input  file  -  line  4.  Verify  the 
location  of  the  depth  file  is  the 
same  as  the  input  file. 

Error  130  -  Opening  Refraction  File. 

readrfrc 

Check  refraction  file  name  in  the 
input  file  -  line  6.  Verify  the 
location  of  the  refraction  file  is  the 
same  as  the  input  file. 

Error  145  -  Input  depth  profile  has  more  data 
points  than  allowed.  Check  depth  profile. 
Program  stopped. 

c_in_dep 

The  maximum  number  of  depth 
points  allowed  is  500.  Modify 
depth  input  file  to  contain  only  500 
depth  values. 

Error  165  -  No  sediment  type  selected  for 
Equilibrium  Profile. 

equilprf 

A  Slope/Sediment  Type  was  not  set 
correctly  in  the  input  file  line  8. 

The  value  must  be  inclusive  of  1-10 

Error  170  -  No  Surf. 

surf 

Check  the  heading  toward  the 
beach  in  the  input  file,  line  7  and 
the  spectrum  input  file.  There  may 
be  no  surf  in  the  area. 

1S1 


Error  180  -  Problem  gridding  to  output  file. 
Program  stops. 

prtoutl 

prt_out2 

Check  that  the  input  depth  profile 
extends  to  the  beach  shoreline  and 
that  the  tide  level  -  line  12  is  not 
too  high. 

Error  185  -  Problem  with  wave  height  values. 

new_brk 

Check  the  input  depth  profile.  The 
data  may  need  to  be  smoothed  due 
to  unusual  slopes.  (Hint:  too  many 
negative  slopes.) 

Error  195  -  Significant  wave  height  outside 
surf  zone  less  than  0.5  ft  -  no  further 
calculations. 

s_nosurf 

Check  the  heading  toward  the 
beach  in  the  input  file  -  line  7. 

Error  200  -  Surf  forecasts  are  for  situations 
when  waves  are  more  important  than  winds. 
This  is  not  the  case  for  input  waves  and 
winds.  Forecasts  may  not  be  valid. 

s_coef  f 

Check  the  input  wave  and  wind 
conditions  in  the  input  file  -  line  1 1 
and  line  12. 

Error  205  -  Water  edge  not  found.  Check  tide 
and/or  depths.  Program  stopped. 

s_tide 

The  input  depth  profile  must  extend 
to  the  beach  including  the  addition 
of  a  tide,  if  specified.  There  must 
be  a  depth  at  either  0.0,  an  onshore 
value,  or  an  elevation. 

Error  210  -  Wave  direction  not  toward  the 
beach  -  no  further  calculations. 

rad_st2 

Check  the  heading  toward  the 
beach  in  the  input  file,  line  7  and/or 
the  directional  wave  spectrum  file. 

Error  215  -  Wave  induced  set-up  not 
converging  to  tolerance. 

setup 

The  input  depth  profile  must  be 
smoothed. 

Error  220  -  Wave  induced  Set-up  is  not 
converging.  Ending  program. 

main_wav 

The  input  depth  profile  must  be 
smoothed. 

Appendix  C.  Flowchart  Symbol  Index 


Terminus  (Start,  Return,  End) 


Process  or  Calculation 


Call  to  Subroutine  or  Function 


a 


Connector 


Keyboard  Input 


/ - \ 


Loop  Limitation 


Preparation  Prior  to  a  Process 
Initialization  of  Variables 


Symbol  Connector 


Yes  Yes  Connector 

I 


—No- 


No  Connector 


Input  or  Output 


Appendix  D.  Acronyms 

CNMOC 

Commander,  Naval  Meteorology  and  Oceanography  Command 

CSCI 

Computer  Software  Configuration  Item 

csu 

Computer  Software  Unit 

DWS 

Directional  Wave  Spectrum 

EOF 

End  of  File 

Hz 

Hertz 

LHS 

Left  Hand  Side  of  Energy  Balance  Equation 

m 

Meter 

N 

Newton 

MSI 

Modified  Surf  Index 

NRL 

Naval  Research  Laboratory 

OAML 

Oceanographic  and  Atmospheric  Master  Library 

ONR  . 

Office  of  Naval  Research 

RHS 

Right  Hand  Side  of  Energy  Balance  Equation 

RSM 

Refraction/Shoaling  Matrix 

SPAWAR 

Space  and  Naval  Warfare  Systems  Command 
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